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ABSTRACT 

We directly measure the evolution of the intergalactic Lyman-a effective optical depth, r e g, over 
the redshift range 2 < z < 4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar 
spectra obtained with the ESI and HIRES spectrographs on Keck, and with the MIKE spectrograph 
on Magellan. This represents an improvement over previous analyses of the Lya forest from high- 
resolution spectra in this redshift interval of a factor of two in the size of the data set alone. We 
pay particular attention to robust error estimation and extensively test for systematic effects. We 
find that our estimates of the quasar continuum levels in the Lya forest obtained by spline fitting are 
systematically biased low, with the magnitude of the bias increasing with redshift, but that this bias 
can be accounted for using mock spectra. The mean fractional error (AC/C tIUC ) is < 1% at z = 2, 
4% at z = 3, and 12% at z = 4. Previous measurements of r c g at z > 3 based on directly fitting 
the quasar continua in the Lya forest, which have generally neglected this effect, are therefore likely 
biased low. We provide estimates of the level of absorption arising from metals in the Lya forest based 
on both direct and statistical metal removal results in the literature, finding that this contribution is 
«6- 9% at z = 3 and decreases monotonically with redshift. The high precision of our measurement, 
attaining 3% in redshift bins of width Az =0.2 around z = 3, indicates significant departures from the 
best-fit power-law redshift evolution (r e g — 0.0018(1 + z) 3 92 , when metals are left in), particularly 
near z = 3.2. The observed downward departure is statistically consistent with a similar feature 
detected in a precision statistical measurement using Sloan Digital Sky Survey spectra by Bernardi 
and coworkers using an independent approach. 

Subject headings: Cosmology: observations, theory — methods: data analysis, statistical, numerical 
— quasars: absorption lines 



1. INTRODUCTION 

The evolution of the intergalactic medium (IGM) 
as traced by the Lyman-a (Lya) forest provides a 
powerful record of the thermal and radiative history 
of the Universe. This power owes to our ability to 
measure the Lya opacity of the IGM as a function 
of redshift, from z = to z > 6, with h igh precision 
from quasar spectra (iPress et alj 119931: iRauch et alJ 



mation within cold d ark matter mode l s (jCen et al.lll99~4l: 
Zhang et all 119951: iHernauist et all 119961: iKatz et all 



19961: iMiralda-Escude et alj |1996t iTheuns et alj 119981 : 
Dave et al.l Il999t ). have been remarkably successful at 



19971 ISongaila et al.lll999h iMcDonald fc Miralda-Escud 
200ll iMeiksin fc Whitef 120031: iTvtler et alJ l2004al 
Songailal 12 004: Bo lton et all [20051 [Kirkman et al l [2005] 
Jena et al.1 120051: iFan et alj l2006ri iBecker et al.l 120071 
Bolton fc Haehnerq l2007h . as well as to the relatively 
simple physics of the Lya forest. In fact, fully and 
pseudo-hydrodynamical cosmological simulations, in 
which the forest arises from absorption by smooth den- 
sity fluctuations imposed on the warm photoionized IGM 
as a natural consequence of hierarchical structure for- 
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4 Department of Astronomy and Astrophysics, UCO/Lick Ob- 
servatory; University of California, 1156 High Street, Santa Cruz, 
CA, 95064, USA. 

5 Jefferson Physical Laboratory, Harvard University, Cambridge, 
MA, 02138, USA. 



reproducing the properties of the absorption observed 
in hi gh-resolution, high signal-to- noise qua s ar spectra 
(e.g-lLu et alJll996tlKirkman fc Tvtlerll 19971 : iKim et all 
I2002air This synergy between theory and observations 
make the Lya forest a particularly compelling probe of 
the diffuse Universe. 

A number of observational results suggest changes in 
the thermal properties of the IGM over the redshift range 
probed by the Lya forest, in particular near z = 3. These 
results fall into four classes: the mean Lya opacity, the 
widths of the Lya absorption lines, the strength of ab- 
sorption by metallic transiti ons, and ultra- v iolet o bser- 
vations of the Hell forest. iBernardi et all (|2003f ) first 
claimed a statistical detection of a feature at z w 3.2 
in the evolution of the Lya forest effective optical depth 
from an analysis of 1061 quasar spectra from the Sloan 
Digital Sky Survey (SDSS; York et al. 2000). These au- 
thors interpreted the observed reduction in the Lya ef- 
fective optical dept h near z = 3.2 as ev i dence for the 
reioni zatio n of Hell (|Theuns et al.ll2002a|) . iRicotti et al.l 
(|2000Q and lSchave et all (|2000f ) both found evidence from 
the Doppler widths of Lya absorption lines for an up- 
ward jump in the temperature of the IGM and a change 
of its equation of s tate t o nea rly isothermal n ear z = 3. 
ISongaila fc Cowid (|1996l ) and ISongailal (|1998h have ar- 
gued for a sudden hardening of the ionizing background 
near z — 3, consistent with the IGM becoming opti- 
cally thin to Hell ionizing photons, based on the ob- 
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served evolution of the CIV/SilV metal-line ratios. Us- 
ing information provid ed by additional metal systems, 
lAgafonova et alj (|2007h argued that the IGM is optically 
thin to Hell ionizing photons at z < 1.8 and measured an 
increasing Hell Lya effect i ve op tical depth from z = 2.4 
to z = 2.9. IVladilo et al.l ()2003l ) used the abundance of 
Arl, which is particularly sensitive to hard photons, in 
damped Lya systems to also argue for a hardening of the 
ionizing background around z = 3. Finally, patchy Hell 
absorption suggestive of the end of Hell reionization has 
been observed in ultra-violet spectra of quasar s HE 2347- 
4342 (jReimers et aBl997llSmette et alJ[2Qpl z = 2.885) 
and Q0320-003 tiJakobsen et al.lll994llHogan et al.l ll997t 
iHeap et al.ll2000L z = 3.286V iFechner et all (|20Q6h inter- 
pret the Hell forest observed in a FUSE spectrum of HS 
1700+6416 (z — 2.72) as arising from residual Hell in its 
post-reionization epoch, while HST/STIS observations of 
QSO 1157+3143 (z = 3) suggest that Hell reionization 
may not be complet e between z = 2.77 and z = 2.97 
(jReimers et al.ll2005t ). 

Each of these observationa l results, however, has 
caveats. iBernardi et al.l (|2003l ) could not estimate the 
Lya optical depth from individual spectra, owing to 
the low resolution and signal-to-noise of the SDSS 
spectra, and instead employed a statistical method. 
While this approach is in principle perfectly valid, 
the procedure is complex and it is difficult to en- 
sure that it is f ree o f systematic effects. Moreover, 
iMcDonald et al.l (j2005j ) also estimated the evolution of 
the mean transmission in the Lya forest from SDSS 
data, usi ng different me t hodol ogies, and did not con- 
firm the IBernardi et al.l (|2003l ) feature, although they 
caution that t heir error es t imate s may not be ac- 
curate . The iRicotti et al.l (|2000f ) and ISchave et al.1 
(2000) temperature and equation of state measurements 
have large error bars and the z ~ 3 temperature 
jump they find is not corrob orated by the ana l yses o f 
IMcDonald et all (l2001h a nd IZaldarriaga et all (hoOlD 
Both the IMcDonald etaH (|200lD and IZaldarriaga et all 
(|2001h results are consistent with an isothermal equa- 
tion of state at z ~ 3, though t he ev idence is weak 
Contr a ry to iSongaila k Cowid (119961) and ISongailal 
(l!998h. iKim et al.l (l2002br ). lBoksenberg et al.l ([2003h and 
lAguirre et al.l (|2004f ) have not found evidence for an 
abrupt change in the CIV/SilV ratio with redshift. The 
z ~ 3 observations of the Hell forest are too scarce to 
draw definitive conclusions on Hell reionization, espe- 
cially if it is inhomogeneous. Moreover, absorption of 
Hell Lya photons is exponentially sensitive to the Hell 
abundance and thus saturates at relatively low densities, 
making it impossible to use the Hell Lya forest to probe 
the heart of Hell reionization. 

On the theoretical side, the rise of quasar l uminos- 
ity function at z ~ 3 (e.g.. iHopkins et aLll2007h is gen- 
erally thojjghtM»_go lia^ 

tion dSokasian et all 12001 iWvithe k Loebl [2003) , since 
stars are expected to emit too few photons hard enough 
(54.4 eV) to ionize Hell. In comparison to the epoch 
of HI reionization, which has recently been the fo- 
cus of detailed theoretical work aimed at predicting 
its timing, ex tent, morphology, and o bservational sig- 
natures (e.g.. [Zaldarriaga et al.l l2004t iFurlanetto et al 



tie attention (for pioneering work, see jCroft et alJll997t 
ISokasian etal]l2002t iTheuns et alll2002aj ) Yet, whereas 
it is unclear whether we have begun to probe the tail of 
HI reionization (for a review of the observatio nal con- 
straints on HI reionization, see iFan e t al. 2006a), or will 
be able to do so in the near future, the observational lines 
of evidence outlined above make a tantalizing case that 
we may be able to infer the details of Hell right now. As a 
phase transition involving 25% of the intergalactic bary- 
onic mass and, possibly, the result of quasar feedback on 
the IGM, Hell reionization is obviously interesting in its 
own right. However, it is also an excellent laboratory to 
test the theoretical framework that has been developed to 
model HI reionization. Indeed, much of the same physics 
underlies the reionization of both HI and Hell. In partic- 
ular, it should be relatively simple to adapt the radiative 
transfer codes developed to simulate HI reionization to 
treat Hell. 

Two lines of attack must be followed to understand 
Hell reionization. First, the caveats on the observational 
results make it clear that further observational work is 
required to clarify the evolution of the thermal and ion- 
ization state of the IGM over the redshift range most 
likely to be affected by Hell reionization, 2 < z < 4. Sec- 
ond, in the absence of large space-based telescopes with 
ultra-violet spectroscopic capabilities to survey the Hell 
forest, the most promising avenue for probing Hell reion- 
ization appears to be in its indirect signatures carried by 
intergalactic hydrogen and metals. In order to draw ro- 
bust and instructive inferences about Hell reionization, 
theoretical work is thus required to quantitatively calcu- 
late these potential signatures. 

Of course, regardless of whether Hell reionization 
is occurring, it is still of great interest to study the 
evolution of the IGM at 2 < z < 4. The mean 
transmission of the Lya forest, for example, probably 
provides the best constraints on t he hydrogen ioniz 
ing background and its evolution 6 dRauch et all 
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Songaila et alj|1999t IMcDonald k Miralda-EscudelL 
Meiksin k Wh ite 2003; Tytler et alJl2004aHBolton et al I 
2005HKirkman et alJl2005b I Jena et al.ll2005h . Indeed, the 
Lya optical depth at any point in the IGM is inversely 
proportional to the photoionization rate. As the ionizing 
background is an integral over all cosmic sources of radi- 
ation, it provides a crucial tracer of star formation and 
quasar activity in the Universe. The thermal state of the 
IGM at these redshifts should also bear relic signatures 
of HI reionization and provide constraints o n that epoch 
(|Theuns et alj|2002bt iHui k Haimanl I2003D . Moreover, 
the mean transmission (F) of the IGM is a very impor- 
tant quantity in the determination of cosmological pa- 
rameters, in pa rticular the matte r power spectrum, from 
the L va forest (ICroft et al.lll998l l2002t IMcDonald et all 



l2^5llLldz~et al.ll2006fl . ISeliak et all (|2003h showed, for 
example, that increasing (F)(z) by 5% at the mean red- 
shift 2.72 changes the best-fit power spectrum amplitude 
by almost a factor of 2, and its slope by 0.2, at the pivot 
point k p = 0.03 s km -1 (see, however. [McDonald et all 



20041 iTrac k Cenll2006HZahn et al.ll2007tlMcQuinn et al 
2007allbl ). Hell reionization has received relatively lit- 



6 In Fauchcr-Gigucrc ct al. (2007), we showed that the existing 
measurements using the proximity effect are likely to biased high 
owing to quasars residing in overdense regions of the Universe. Ob- 
servational difficulties, such as uncertainties on quasar systematic 
redshifts, and temporal variability of quasar luminosities, also make 
it challenging to obtain robust proximity effect measurements. 
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120051 who bypass the use of an external constraint on (F) 
in their estimation of the linear theory power spectrum 
from SDSS data). 

In this paper, we make an important step in the obser- 
vational direction by presenting a direct precision mea- 
surement of the effective Lya optical depth and its evo- 
lution over the redshift range 2 < z < 4.2 from a large 
sample of 86 high-resolution, high signal-to-noise quasar 
spectra obtained obtained with the ESI and HIRES spec- 
trographs on Keck and with MIKE on Magellan. Our 
focus here is on the presentation of model-independent 
observational results with robust uncertainty estimates. 
As such, we provide detailed numerical tables and fitting 
formulae, and leave theoretical interpretations to sepa- 
rate work. 

We describe our data set in SJSJ Our measurement of 
the Lya effective optical depth and its statistical error 
analysis are presented in Sj3j We also present estimates 
of the bias arising from the difficulty of estimating quasar 
continua at high redshifts and of the level of absorption 
from metals in the Lya forest in this section. The sys- 
tematic uncertainty on the continuum and metal correc- 
tions is estimated in 21 In $5j we fit simple functions 
to our measurement, corrected for metal absorption and 
not, and compare our results with previous work. We 
conclude in Sj6] by outlining the physical effects that may 
give rise to the feature we find evidence for near z = 3.2 
and their degeneracies, leaving theoretical interpretation 
for future work. A detailed appendix describes a series 
of tests for systematic errors. 

Throughout this paper, we assume a cosmology with 
(Sim, Sl b , Sl A , h, a 8 ) = (0.27, 0.046, 0.73, 0.7, 0.8), 
close to the parameters inferred from the Wilkinson Mi- 
crow ave Anisotropy Probe three- year data (jSpergel et alJ 
12001 . 

2. DATA SET 

Our data set consists of quasar spectra ob- 
tained with the High Resolution Echelle Spectrometer 
( Vogt et aLlll99 4, HIRES) and Spectrograph and Imager 
(jSheinis et al.l l2002. ESI) spectrographs on the Keck tele- 
scopes and with the M agellan Inamori Kyocera Echelle 
(jBernstein et al.l 120031 MIKE) spectrograph on Magel- 
lan. 

For nearly all observations, the HIRES spectra were 
acquired using either a 0.8" or 1.1" wide decker 
(FWHM ps 6 and 8 km s" 1 , respectively) and the ESI 
observations were carried out with the 0.5" or 0.75" slit 
(FWHM ps 33 and 44 km s _1 , respectively). In gen- 
eral, the signal-to-noise ratio (S/N) of these spectra is 
> 15 per pixel (2 km s _1 pix -1 for HIRES and 11 km 
s- 1 pix- 1 for ESI). Most of the ESI and HIRES spectra 
were originally obtained to study damped Lya absorbers 
(DLAs). A subset of these has been public ly released 
to the community by IProchaska et al.l (l2007f) a n d wer e 
previously summarize d in iProc haska fc Wolfd (119990 . 
IProchaska et all (|2001[ ). and IProchaska et al.l (I2003D. A 
few ot her ESI spectra were analyzed in lO'Meara et all 
(2007). This sample is supplemented by a number of new 
ESI spectra also selected for DLA studies (Prochaska et 
al., in prep.). Seven other HIRES spectra were provid ed 
by M. Rauch and W. Sargent (e.g. lRauch et al.1119971 ). 

The MIKE spectrograph has a dichroic optical ele- 
ment splitting the beam into blue and red arms, with 
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Fig. 1. — Examples of quasar spectra used in this study. The top 
panel shows a spectrum obtained with Keck/HIRES, the middle 
panel a spectrum obtained with Keck/ESI, and the bottom panel 
one obtained with Magellan/MIKE. The spectra were normalized 
by d ividing them by the estimated continuum level (_F = 1) (see 
i]3.4l l. indicated by the thin horizontal line. The ESI spectrum has 
lower resolution and fewer pixels. 

FWHM « 11 km s" 1 and 14 km s _1 for the blue and 
red sides, respectively, for a 1.0" slit. Our MIKE spectra 
have S/N > 10 per pixel and were analyzed fo r supe r 
Lyman-limit systems (SLLS) in lO'Meara et~aT1 (|2007h . 
Figure [1] shows examples of quasar spectra obtained with 
each instrument used in this study. Tables [TJ ^ and [3] 
lists the quasars in the HIRES (16), ESI (44), and MIKE 
(26) data sets, respectively, with their redshifts. Figure 
[2] illustrates the quasar redshift distributions, as well as 
the corresponding path lengths of Lya absorption, in his- 
togram form. 

3. MEASUREMENT OF THE EFFECTIVE Lya OPTICAL 
DEPTH 

We now measure the Lya opacity of the IGM as a func- 
tion of redshift from the data set described in $2] We be- 
gin by defining the quantities of interest and describing 
our methodology, and present the measurement itself, in 
^3.11 Masks which are applied to the data are described 
in i)3.2l We discuss the calculation of error bars in ^3.31 
In §3.4[ we test for systematic effects in estimating the 
continuum level of the quasar spectra, and calculate the 
redshift-dependent bias arising from the increasing level 
of absorption with increasing redshift, which we use to 
correct our measurement. We estimate the level of ab- 
sorption in the Lya forest arising from metals in i j3.5l 

3.1. Definitions and Methodology 

Let i 7 'abs(A) be a quasar's absolute flux and C(A) be its 
continuum (unabsorbed) level as a function of observed 
wavelength A. The corresponding redshift of Lya absorp- 

^Lya 



tion is just z LyQ , = A/A Lya - I . when- A T = 1216 A. 
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Fig. 2. — Histograms showing the redshift distribution of the 
quasars (top panel) used in this study and the corresponding proper 
path length distribution of Lya absorption (bottom panel). In each 
panel, the red histogram corresponds to HIRES data only, the blue 
one to HIRES+ESI data, and the black one to the complete data 
set (HIRES+ESI+MIKE). 



TABLE 1 
HIRES Quasar 
Sample 



Name z em 



Q0000-26 a 


4.11 


BR0019-15 


4.53 


Q0336-01 


3.20 


Q1107+49 a 


3.00 


Q1209+09 


3.30 


Q1210+17 


2.54 


Q1215+33 


2.61 


Q1422+23 a 


3.62 


Q1425+60 


3.20 


Q1442+29 a 


2.67 


Q1759+75 


3.05 


Q2206-19 


2.56 


Q2237-06 a 


4.55 


Q2231-002 


3.02 


Q2343+12 a 


2.52 


Q2359-02 


2.80 



a Provided by 
M. Rauch and 
W. Sargent (e.g ., 
IRauch et. aT1IT59"7l) . 



Dropping the Lya subscript on z, we define the trans- 



TABLE 2 
ESI Quasar Sample 



Name z em 



PSS0007+2417 4.05 

SDSS0013+1358 3.58 

PX0034+16 4.29 

Q0112-30 2.99 

SDSS0127-00 4.06 

PSS0134+3307 4.52 

SDSS0139-0824 3.02 

SDSS0142+0023 3.39 

PSS0209+0517 4.17 

SDSS0225+0054 2.97 

BRJ0426-2202 4.32 

PSS0808+52 4.45 

SDSS0814+5029 3.88 

Q0821+31 2.61 

SDSS0844+5153 3.20 

SDSS0912+5621 3.00 

SDSS0912-0047 2.86 

Q0930+28 3.42 

BQ1021+30 3.12 

CTQ460 3.13 

HS1132+22 2.89 

SDSS1155+0530 3.48 

PSS1159+13 4.07 

PSS1248+31 4.35 

PSS1253-02 4.01 

SDSS1410+5111 3.21 

PSS1432+39 4.28 

HS1437+30 2.99 

PSS1506+52 4.18 



PSS1535+2943 3.99 
SDSS1610+4724 3.22 



PSS1715+3809 4.52 

PSS1723+2243 4.52 

PSS1802+5616 4.18 

FJ2129+00 2.96 

PSS2155+1358 4.26 

Q2223+20 3.56 



SDSS2238+0016 3.47 

PSS2241+1352 4.44 

SDSS2315+1456 3.39 

PSS2323+2758 4.18 

FJ2334-09 3.33 

PSS2344+0342 4.24 

SDSS2350-00 3.01 

mission, or norm alized flux, a3 

Let (F) be the ensemble average over quasars and further 
define the effective optical depth as 

reff(*) = - In [<*">(*)] (2) 

(e.g-. lRauch et alJ H.997). This quantity, being a function 
of (F) only, has the convenient property of being inde- 
pendent of noise, provided the noise fluctuations have a 
symmetric distribution about zero and sufficiently many 
independent pixels are averaged over. Both conditions 
are satisfied for our large sample. Moreover, it is stan- 
dard in the literature, which will allow us to directly 
compare our results with published measurements. 
The measurement consists of the following steps: 

1. Estimate the continuum level in the Lya forest 
of each quasar spectrum and calculate the corre- 
sponding F; 

2. apply masks to the data: quasar proximity regions, 
higher-order Lyman series forests, SLLS, DLAs, 
and associated metal lines; 
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TABLE 3 
MIKE Quasar Sample 



Name z em 



nmm 304 


3 


14 




2 


. ' ) (.) 


UFO'340 9fil 9 


Q 
O 


OS 


SDSS091 9-1-0^47 




.25 


Fl ^1 0Q4 9 J_ 04 9 9 


■? 

o 


97 
. Z ( 






1 n 


Cjrmm 09^4-04^9 




24 


ni 1 no 9fSzL 


2 


14 








OI 994 flSI 9 
V^lZZ^t-UolZ 


o 
Z. 


"1 -1 


cin^im 94Q m 


■3 

O. 


AH) 


^Fl^Cll ^Q-LO^S 
oLJijo looy-rUO^o 


o 

Z. 


07 

. y * 


HE1347-2457 


2 


.60 


Q1358+1154 


2 


.58 


SDSS1402+0146 


4 


.19 


SDSS1429-0145 


3 


.42 


Q1456-1938 


3 


.16 


SDSS1621-0042 


3 


.70 


PKS2000-330 


3 


.78 


Q2044-1650 


1 


.94 


Q2126-158 


3 


.28 


HE2156-4020 


2 


.53 


HE2215-6206 


3 


.32 


HE2314-3405 


2 


.94 


HE2348-1444 


2 


.93 


HE2355-5457 


2 


.93 



3. bin the remaining Lya forest pixels in redshift in- 
tervals; 

4. for each redshift bin, average the pixels to estimate 
(F) and calculate r c g using Equation @; 

5. correct this "raw" measurement for redshift- 
dcpcndcnt continuum estimation error and, option- 
ally, for metal absorption. 

The continuum level of each spectrum is estimated by 
fitting a cubic spline through its transmission peaks for 
the ESI and MIKE spectra. To do so, we use the IDL 
graphical user interface x_continmim developed by J. X. 
Prochaska. More precisely, we manually identify peaks 
of transmission near unity by looking for plateaux con- 
sisting of several consecutive pixels with approximately 
the same transmission level and force the spline to pass 
through these points. In practice, these plateaux have a 
finite thickness owing to the noise level and so we place 
the points at the centers of the noise bands. Additional 
spline points are then placed and adjusted to approxi- 
mate the most featureless continuum consistent with the 
identified transmission peaks. In general, the number of 
transmission peaks that we are able to identify decreases 
with redshift. The robustness and accuracy of this man- 
ual spectrum-by-spectrum continuum estimation method 
are tested in £13.41 and Appendix [Bl For the HIRES data, 
the continua were similarly estimated as part of the pro- 
cedure to join the different echelle orders. Figures [1] and 
1181 show examples of continuum fits on real and mock 
spectra, respectively. The masks, detailed in the next 
section, are necessary to ensure that the measurement 
is uncontaminated. Pixels, which actually have a finite 
width, are labeled by their central redshift in perform- 
ing the binning. Any given pixel therefore falls in the 
bin which contains > 50% of it. Errors are calculated 
as in £13.31 Our measurement from the complete data 
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set (HIRES+ESI+MIKE) in Az = 0.2 redshift bins, raw 
and corrected for continuum bias as in H3.41 is shown in 
Figure [3J The numerical values are provided in Table 21 
in which we compile all of our r e g measurements, cor- 
rected for metal absorption (see §3.5() and not, both in 
Az = 0.1 and Az = 0.2 redshift bins. 

3.2. Masks 
3.2.1. Definition of the Lya Forest 

To avoid contamination from absorption by higher 
Ly man-series resonances, we consider only those pixels 
between the rest- frame wavelengths of Ly/3 (1025 A) and 
Lya (1216 A). Further, we mask the "proximity region" 
of each quasar, in which the absorption is biased low with 
respect to the cosmic mean owing to the ionizing radi- 
ation of the quasar itself. Specifically, we mask the 25 
proper Mpc nearest to each quasar. The Hi-ionizing radi- 
ation of bright quasars typically equals the cosmic mean 
at dis tances ~ 5 proper Mpc fe.g. jFaucher-Giguere et ail 
2007), so that our masks are very conservative and little 
contamination should remain. 

3.2.2. Bad and Low-S/N Pixels 

Some quasar spectra contain bad pixels, resulting in 
gaps in the data. These bad pixels are identified by 
carefully inspecting each spectrum and are then masked. 
The sensitivity of the spectrographs generally degrades 
at shorter wavelengths. The correspondingly low S/N 
makes reliable continuum estimation challenging. We 
thus keep only the portions of the spectra where we 
judge that the continuum can be reliably estimated, cor- 
responding to S/N > 10 per pixel. 

3.2.3. Damped Absorbers 

Most of the quasar spectra used here were obtained 
as part of studies of SLLS and DLAs. Since these sys- 
tems individually contribute large equivalent widths to 
the observed absorption, we mask these systems to en- 
sure that our measurement of the Lya forest mean trans- 
mission is unbiased. Specifically, we mask systems with 
log 10 Nm ^ 19-0 and associated metal lines. We show in 
Appendix IA.3I that these systems contribute only about 
2% of the total absorption in the Lya forest in an unbi- 
ased sample. This is the amount by which the absorption 
is underestimated by masking these systems in our biased 
sample. The width of each mask is chosen such that all 
data around the SLLS or DLAs where the transmission 
in the presence of that system only would be < 99% is 
ignored. Beyond masking the absorption from the SLLS 
or DLAs, this helps to keep the error on the estimated 
continuum owing to the damping wings < 1%. In Ap- 
pendix IA.31 we increase the width of the masks on the 
SLLS and DLAs by a factor of two and conclude that our 
measurement is robust to the choice of width. For the 
metal lines, we mask the most common neutral species 
(Mgl, OI, etc.), low-ionization ions - which are dominant 
in HI absorbers - (Fell, Sill, CII, OI, NI, A1II, etc.), in- 
termediate ions (A1III, Felll, CHI, Silll, etc.), and high- 
ionization ions (SilV, CIV, NV, OVI, etc.). We place a 
mask of width 500 km s _1 centered on each metal line as- 
sociated with a SLLS or DLA in our sample. The width 
of metal lines is generally dominated by Doppler broad- 
ening owing to peculiar motion, for which the median is 
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Fig. 3. — Measurement of T e g vs. 2 in Az = 0.2 bins from the complete data set. The red circles sho w the raw measurement obtained 
from the continuum fits. The bl ue s quares show the measurement after correction for continuum bias as in t|3,4l Error bars are one standard 
deviation and calculated as in 33.31 The continuum-corrected data points are centered on the same redshifts as the raw measurement, but 
have been slightly offset to the right for graphical clarity. Numerical values are provided in Table [4] 

verified that the error estimates have converged for our 
choice of segment length. In the limit of an infinite 
number of segments, the central limit theorem ensures 
that the uncertainty on (F) is Gaussian. Figure 0] shows 
histograms of (F)i for the redshift bins considered. 
The histograms are close to Gaussian, so that the 
convergence to Gaussianity of the error on (F) should 
be rapid, and may thus confidently be expected to have 
been attained. The histograms of Figure 0] also show 
no sign of outliers, arguing that our measurement is 
unaffected by them. 



-80 km s" 1 (jProchaska k. Wolfd fl99l . Our masks are 
thus sufficiently large to completely masks these lines 
and allow for an uncertainty — 200 km s _1 in each direc- 
tion on the redshift of the lines. We show in Appendix 
IA.3I that masking metal lines associated with DLAs has 
practically no effect on our measurement. 

3.3. Statistical Error Bars 

To estimate the statistical uncertainty on our measure- 
ment of (F), we proceed as follows. For any given redshift 
bin, we first concatenate the data from the different spec- 
tra. We then break the total data stream in segments of 
length 3 proper Mpc, where the length is chosen such that 
the mean transmission calculated from t he i-th segment, 
(F)j, is ind ependent of its n e ighbo rs. iMcDonald et al.1 
( 2000) and IMcDonald et all (|2006l) measured the flux 
correlation function and power spectrum in the Lya for- 
est an d found small correlat i ons o n scales >1 proper 
Mpc. iFaucher-Giguere et al.l (|2007f l computed the op- 
tical depth correlation coefficient p from a simulation of 
the Lya forest and found p < 1% for r > 3 proper Mpc 
over the redshift range 2 < z < 4. The error on (F) is 
then estimated as 



°{F) = 



N 



(3) 



where ct/f\ ( is the standard deviation among the in- 
dividual (F)ii and N is the number of independent 
segments. This estimate simply makes use of the fact 
that (F) is the mean of the (F)i, which are mutually 
independent, and of the usual result that the error on 
the mean of independent measurements is just equal 
to their standard deviation divided by y/~N. We have 



In 21 we estimate the systematic error budget, in- 
cluding contributions from the continuum and metal 
corrections (see below). 

3.4. Continuum Bias 

Our measurement of (F) is a direct function of our 
continuum estimates, and we must therefore ask how ac- 
curate these are. 

Random continuum errors are included in the statis- 
tical error estimates described We thus need only 
further consider potential systematic biases: finite res- 
olution and S/N of the data, and a possible redshift- 
dependent bias arising from increasing Lya absorption 
with increasing redshift. 

Finite resolution and S/N both make accurate con- 
tinuum fitting more challenging. In the case of finite 
S/N, the pixel-wise "peaks" in the observed spectra in 
general do not exactly correspond to a transmission of 
unity, owing to the noise contribution. Since the noise 
fluctuations are symmetric about zero and we attempt to 
identify the transmission peaks with the centers of noise 
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TABLE 4 

Measurement of r e g, with and without continuum and metal corrections, and in redshift bins of width Az = 0.2 and 

Az = 0.1 



z R,aw a Cont. corr. b Schaye metal corr. c Kirkman metal corr. c o> cff ,stat d <JY cff ,cont c °V,, ff , metal Q"T cff .tot s 

Az = 0.2 — — — 



2.0 


0.138 


0.146 


0.127 


0.115 


0.018 


0.002 


0.012 


0.023 


2.2 


0.176 


0.187 


0.164 


0.157 


0.013 


0.003 


0.014 


0.020 


2.1 


0.214 


0.229 


0.203 


0.200 


0.009 


0.004 


0.016 


0.019 


2.6 


0.258 


0.280 


0.251 


0.251 


0.010 


0.006 


0.018 


0.022 


2.S 


0.330 


0.360 


0.325 


0.332 


0.012 


0.008 


0.021 


0.026 


.3.0 


0.384 


0.423 


0.386 


0.396 


0.014 


0.010 


0.022 


0.029 


3.2 


0.399 


0.452 


0.415 


0.425 


0.017 


0.012 


0.022 


0.031 


3.4 


0.546 


0.615 


0.570 


0.589 


0.019 


0.017 


0.027 


0.039 


.3.6 


0.677 


0.766 


0.716 


0.742 


0.030 


0.023 


0.030 


0.051 


3.8 


0.770 


0.885 


0.832 


0.861 


0.025 


0.028 


0.031 


0.051 


4.0 


0.839 


0.986 


0.934 


0.963 


0.037 


0.034 


0.031 


0.064 


4.2 


0.926 


1.113 


1.061 


1.090 


0.091 


0.041 


0.031 


0.121 



Az = 0.1 



2.0 


0.148 


0.156 


0.135 


0.125 


0.022 


0.002 


0.012 


0.027 


2.1 


0.130 


0.139 


0.121 


0.109 


0.017 


0.002 


0.011 


0.022 


2.2 


0.173 


0.184 


0.161 


0.154 


0.016 


0.003 


0.013 


0.023 


2.3 


0.230 


0.243 


0.215 


0.213 


0.018 


0.004 


0.017 


0.026 


2.4 


0.184 


0.199 


0.177 


0.170 


0.010 


0.004 


0.014 


0.018 


2.5 


0.238 


0.257 


0.229 


0.228 


0.012 


0.005 


0.017 


0.022 


2.6 


0.269 


0.291 


0.260 


0.262 


0.013 


0.006 


0.018 


0.024 


2.7 


0.284 


0.309 


0.278 


0.281 


0.015 


0.006 


0.019 


0.026 


2.8 


0.334 


0.363 


0.328 


0.335 


0.018 


0.008 


0.021 


0.030 


2.9 


0.339 


0.373 


0.339 


0.346 


0.016 


0.008 


0.021 


0.029 


3.0 


0.401 


0.441 


0.402 


0.413 


0.021 


0.010 


0.023 


0.035 


3.1 


0.405 


0.450 


0.412 


0.424 


0.023 


0.011 


0.023 


0.036 


3.2 


0.402 


0.454 


0.418 


0.428 


0.023 


0.012 


0.022 


0.036 


3.3 


0.419 


0.479 


0.442 


0.453 


0.025 


0.013 


0.022 


0.038 


3.4 


0.533 


0.602 


0.558 


0.576 


0.025 


0.017 


0.026 


0.043 


3.5 


0.685 


0.763 


0.710 


0.738 


0.029 


0.022 


0.032 


0.051 


3.6 


0.666 


0.755 


0.705 


0.731 


0.046 


0.023 


0.030 


0.065 


3.7 


0.718 


0.820 


0.768 


0.795 


0.042 


0.025 


0.030 


0.062 


.3.8 


0.753 


0.868 


0.817 


0.844 


0.033 


0.028 


0.031 


0.057 


3.9 


0.802 


0.931 


0.879 


0.908 


0.029 


0.031 


0.031 


0.056 


1.0 


0.782 


0.929 


0.880 


0.906 


0.058 


0.032 


0.029 


0.082 


1.1 


0.953 


1.118 


1.063 


1.095 


0.082 


0.040 


0.033 


0.110 


4.2 


0.893 


1.080 


1.030 


1.058 


0.101 


0.040 


0.030 


0.132 



a Estimated directly from the continuum fits as in i]3.1l b Raw measurement corrected for the redshift-dependent bias in the continuum fits 
as in H3,4I C Same as the measurement corrected for continuum bias, but with the contribution of metals to absorption in the Lyc* forest 
corrected as in i]3.5l d Statistical error, estimated as in i|3.3l c Systematic error arising from the continuum correction, estimated as in 
5 14.11* Systematic error arising from the metal correction, estimated as in i|4.2l g Quadrature sum of the statistical and systematic errors 

(S3)- 

bands, finite S/N is not expected to contribute a signif- 
icant systematic bias, provided S/N > 10 per pixel. To 
verify this, we plot the ratios of (F)(z) as estimated from 
every particular quasar intersecting a given redshift bin, 
(F)(z)qso, to (F)(z) for the entire bin (as calculated 
from all the quasars contributing to it) versus the mean 
noise in the particular quasar's data in the bin in Figure 
[H We use redshift bins of width Az = 0.2. These ratios 
should average to unity in the absence of systematics. No 
significant trend with noise is seen, thus arguing against 
noise-dependent systematic effects. 

In the case of finite resolution, the peaks are smoothed 
and their height may therefore be underestimated. Is 
this effect important for our data? The pressure of the 
baryons smoothes their spatial distribution on a scale of 
order the Jeans scale, 



' TT"fkT 

Gpfirrip 



(4) 



where 7 is the adiabatic index of the gas, T is its temper- 
ature, p is the total mass density, p the mean molecular 
weight of the gas, and m p is the proton mass. To convert 



to km s , we multiply A, 7 by H(z) = ^/SitCp/Z (the 
Friedmann equation for a flat universe). For 7 = 1.62 
(valid in the limit of early reionization; iHui &: Gnedinl 
Il997f ). p = 0.59 (for a fully ionized gas consisting of 75% 
hydrogen and 25% helium by mass) and T = 20, 000 K 
(iRicotti et all 120001: ISchave et all 120001 : iMcDonald et all 
I200H IZaldarriaga et al.ll2001f ). A 7 « 103 km s" 1 , larger 
than the resolution of each spectrograph used in this 
study (c.f. Ej2j) . If the transmission peaks correspond to 
density troughs of size at least the Jeans scale, this sug- 
gests that all of our spectra will resolve them, and hence 
that our measurement will not be affected by any signifi- 
cant resolution-dependent bias. This argument, however, 
has caveats. 

The correct baryonic smoothing scale is in fact not 
equal to the Jeans scale, but instead is expected to 
be smaller, owing to the finite time necessary for the 
baryons to respond to the heat injected into the IGM 
durin g HI reionization at redshift > 6 (Gn edin fc Huil 
1998). On the other hand, thermal broadening redshift- 
space distortions will additionally smooth the spectra 
on a scale comparable to the Jeans scale (for details 
about redshift-space distortions in the Lya forest, see e.g. 
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Fig. 4. — Histograms of the segment- averaged transmissions (F)i. Each panel corresponds to a redshift bin of width Az = 0.2. The text 
in the upper left corners indicates the total number of 3 proper Mpc segments in the bin, the mean (-F);, and their standard deviation. 
The 1<t uncertainty on (F) = {{F)i) is estimated as a if) = CT (F).J y J r N ( ^3.31 1. The solid curves shows the Gaussian function with same 
mean and standard deviation as the histograms. 

ally below the cosmic mean, reducing the relevant local 
^ . smoothing scale, and may not be definitively known. To 
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Fig. 5. — Ratios of (F)(z) as estimated from every particular 
quasar intersecting a given redshift bin, (F)(z)qso, to (F)(z) for 
the entire bin (as calculated from all the quasars contributing to 
it) versus the mean noise on F in the particular quasar's data 
in the bin. The crosses show averages over noise bins of width 
0.05. The ler vertical error bars are equal the standard deviation 
of the points in the bin divided by the square root of the number 
of points. No significant trend with noise is seen, arguing against 
noise-dependent systematic effects. 



iFaucher-Giguere et alj [20071. The Lya forest transmit- 
ted flux is also measure d to have non-zero powe r on scales 
smaller than the Jeans (McDon ald et al.|[2~000 f). Further, 
the IGM temperature in underdense regions is gener- 
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Fig. 6. — Mean ratio of {F)(z) as estimated from every particular 
quasar intersecting a given redshift bin, {F)(z)qso> t° (F)(z) for 
the entire bin (as calculated from all the quasars contributing to 
it) versus the instrument with which each spectrum was obtained. 
The FWHM of each spectrograph is indicated on the figure (Sj2j. 
The true resolution varies from spectrum to spectrum, depending 
on the width of the slit used; here we show an optimistic, yet rep- 
resentative, value for each spectrograph. The lcr vertical error bars 
are equal the standard deviation of the points in the bin divided by 
the square root of the number of points. No significant instrument- 
specific bias or trend with resolution is seen, arguing against such 
systematic effects. 

test for resolution- or spectrograph-dependent systemat- 
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ics from the data itself, we compare in Figure [6] the mean 
ratio of (F)(z) as estimated from every particular quasar 
intersecting a given redshift bin, (F}(z)qso, to {F}(z) 
for the entire bin (as calculated from all the quasars con- 
tributing to it) versus the instrument with which each 
spectrum was obtained. Again, no significant bias or 
trend with instrument or resolution is seen. 

Finally, as the redshift increases, the IGM density in- 
creases owing to the cosmological expansion. As a re- 
sult, the Lya absorption also increases with redshift and 
transmission peaks reaching unity become increasingly 
rare, potentially causing us to underestimate the con- 
tinuum level more at higher redshifts. In Appendix [Bj 
we detail how we use mock spectra to quantify the mean 
fractional difference between the true and estimated con- 
tinua, (AC/Ctruo) (where AC = C cs t — Ct mc ), as func- 
tion of redshift. We find that, for hydrogen background 
photoionization rates (T bk s) taken from the literature, 
AC/CUc = 1-58 x 1CT 5 (1 + z) 5 63 to sub-percent sta- 
tistical precision over the redshift range 2 < z < 4.5, 
as illustrated in Figure [7] We note, however, that this 
continuum correction itself is subject to some systematic 
bias, since it depends on assumptions (particularly re- 
garding T bkg ) used to generate the mock spectra. This 
potential systematic bias is explored in detail in the Ap- 
pendix. In the Appendix, we also show that if r^? is 
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Fig. 7. — Mean fractional difference between the estimated and 
true continua as a function of redshift. The points at redshift 
intervals Az = 0.5 were calculat ed b y comparing blind continuum 
fits on mock spectra (see Figure [18} to true continua for 10 mock 
spectra at each redshift, as described in i|3,4l and Appendix|E] The 
small error bars indicate the la errors on the means. The redshift 
trend is well approximated by the best-fit power law AC/Ctrue = 
1.58 X 10 — 5 (1 + z) 5 63 over the redshift range covered. 

the effective optical depth as estimated as above, and 
r true j g "tj ue " va j ue after correcting for the redshift- 
dependent continuum bias, then 



'off 



111 



1 



AC 

Ctrue 



(5) 



We use this expression to correct our measurement for 
continuum error; the continuum-corrected measurement 



is compared to the raw measurement in Figure[3l We will 
henceforth generally consider this continuum-corrected 
measurement. 

Of course, this effect does not affect measurements 
based on extrapolating the continuum from redward 
of the Lya line, as usually done in measurements 
based on lower resolution data for which direct local 
continuum estimation is impossible (e.g., iPress et al.1 
[l^lBernardi et al.ll200l iMcDonald et al.ll2005r) . Such 
extrapolations, however, likely miss local variations 
in the continuum from quasar to quasar so that local 
estimation is advantageous for high-resolution data. 

3.5. Metal Absorption 

Metals also contribute to the total absorption in the 
Lya forest. Consider the transmission at observed wave- 
length A, F(X), in a given spectrum. Generically labeling 
metals by m, 



F(\) = F a xJ[F„ 



it 



(6) 



where r a is the Lya optical depth and r m is the optical 
depth owing to metal to. 7 Because each metal transition 
has a different rest-frame wavelength, absorption mani- 
festing itself at observed wavelength A occurs in physi- 
cally separated regions of the IGM. The F m are therefore 
statistically independent, so 



(F(X)) = (F a ) x IJ<F„ 



From this, 



T e ff = Toff., 



, r eff,; 



(7) 



(8) 



where we have simply generalized the effective optical 
depth definition (eq. [2]). Equation ([8]) states that the 
measured effective optical depth is actually the sum of 
a contribution from pure Lya and contributions from 
metals. 

The easiest quantity to model, either analytically 
or by simulation, is T c ff,a. 

Can we correct our measurement to estimate T e ff iQ ? To 
a certain extent, it is possible to separate metal absorp- 
tion from Lya absorption. We explore two ways, one 
direct and one statistical, based on published measure- 
ments; 

iSchave et al.l (|2003| ) measured r e ff from 19 high-quality 
spectra obtained with the UV- Visual Echcllc Spectro- 
graph (UVES) instrument on the Very Large Telescope 
(VLT) and HIRES over approximately the same redshift 
range as us (2 < z < 4.2). These authors attempted 
to directly identify contaminating metal lines by search- 
ing by eye for absorption features corresponding to the 
redshifts of strong HI lines and to the redshifts of the ab- 
sorption systems redward of the quasar's Lya emission 
line. They fitted a power law of the form 

logr ofr = logTo + alog[(l + z)/4] (9) 

to both the data with all the metals and with metals re- 
moved. When metals are kept, they found (log To, a) = 

7 For the purpose of this discussion, different metal ionization 
states are viewed as different species. 
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(-0.40 ±0.02, 3.40 ±0.23), and when metals are masked 
(log T , a) = (-0.44 ± 0.01, 3.57 ± 0.20). We must cau- 
tion in passing that the la error bars on the lSchave et all 
(2003) r e fi measurements appear to be significantly un- 
derestimated based on comparing the scatter between 
the different data points to their reported errors (see 
their Figure 1 ). As suming that the relation between the 
ISchave et alj (|2003[ ) fits accurately represents the metal 
contribution to the absorption in the Lya forest, we can 
calculate a correction to our more precise T e g measure- 
ment: 



'HD, metals \ / 

(10) 

In Figure [HI we compare our r c g measurement with 
T c g t a inferred using this correction. The relative metal 
correction to r e g decreases with increasing redshift, as ex- 
pected if the IGM is monotonically enriched with time 8 : 
13% at z = 2, 9% at z = 3, and 5% at z = 4. 

An alternative to directly identifying metal lines in the 
Lya forest is to statistically remove them using mea- 
surements of "pure" metal absor p tion re dward of the 
Lya emission line. iTvtler et al.1 (|2004bf ) applied this 
idea in two ways at redshift ~ 1.9. First, they mea- 
sured the absorption in the rest-frame wavelength range 
1225-1500 A in their sample of spectra obtained with 
the Kast spectrograph on the Lick Observatory 3 m tele- 
scope. As more metal lines contribute to the absorp- 
tion at shorter rest-frame wavelengths, they extrapo- 
lated to estimate the metal absorption in t he Lya forest 
wavelen gth range 1070-1170 A. As a check. ITvtler et al.1 
(20041)) also summed the equivalent widths of metal ab- 
sorption l ines for 26 quasars w ith 1.7 < z em < 2.3 tab- 
ulated bv lSargent et al.l <l!988f) to es t imate the amount 
of absorption by metals. ITvtler et al.l (|2004ph considered 
the quantity DA = 1 — (F) instead of r e g, but note that 
(F) = exp (—r ff) « 1 — T c ff and hence T e ff ~ DA for 
T c ff <C 1, which is a good approximation at z ~ 2. Defin- 
ing DM to be the contribu tion to DA owing to met- 
als only. ITvtler et all (|2004bD found DM = 0.158 ± 0.13 
from their Kast analysis (their DM2) an d 0.165 ± 0.22 
(their DM6) using the lines tabulated by I Sargent et al.l 
(1988) at observed wavelength 4158 A. The two val- 
ues thus agree very well. Building upon the work of 
ITvtler et all (l2 004b1. Kirkman et~al1 ([20051 ) included all 
52 lSargent et a l. ( 1988) quasars to estimate the metal ab- 
sorption over the extended redshift range 1.7 < z < 3.54. 
They estimated the metal contribution as a function of 
both rest-frame (A r ) and observed wavelengths (A D ). In 
their notation: 

DM5(X r ) = 0.01564 - 4.646 x 10 _5 (A r - 1360 A) (11) 

and 

DM6(A D ) = 0.01686- 1.798 x 10~ 6 (A o - 4158 A). (12) 

Although the observed wavelength is the relevant quan- 
tity for the redshift of Lya absorption, the trend with 
rest-frame wavelength serves to approximately account 

8 Note, however, that the metal contribution to the absorption 
(i.e., r e ff) need not strictly follow the abundance of the metals, 
since the amount of absorption also depends on the ionization state 
of the metals, which in turn depends on parameters like the inten- 
sity and hardness of the radiation background. 



for the growing collection of metals contributing absorp- 
tion at shorter wavelengths. We first estimate the metal 
contribution a t the mean rest-frame wavelength of the 
ISargent et al.l (|1988l ) sample, (A r ) = 1360 A, and at 
Lya redshift z using Equation (fT2|) : 

T e s,m( z = Ao/1216 A - 1; A r = 1360 A) = 

-ln[l-£»M6(A )]. (13) 

We then use Equation (fTTj) to extrapolate T e g t7n to a 
wavelength approximately in the center of the Lya forest 
region, 1120 A: 

T e g,m(z; K = 1120 A) _ In [1 - DM5(X r = 1120 A)] 
T e ff,m(z; A r = 1360 A) ~ In [1 - DM5(\ r = 1360 A)] 
= 1.723. (14) 

In Figure [HI we compare our T c ff measurement cor- 
rected only for continuum bias with both T c s yCC inferred 
using this statistical method correction and wit h the 
direct method based on the ISchave et al.l (|2003t ) fits. 
The fraction of T ff contributed by metals is 21% at 
2 = 2, 6% at z = 3, and 2% at z — 4 according to 
the statistical metal correction. Our metal-corrected 
Toff, a estimates based on removing metals directly and 
statistically agree well within la at each redshift con- 
sidered, thus suggesting that the metal correction using 
either method is accurate to the level of our statistical 
error bars, or better. A more detailed estimate of the 
uncertainty on T e ff, Q introduced by the metal correction 
is provided in JJU 

Numerical values are tabulated in Table [H 



4. SYSTEMATIC ERRORS 

In the Appendix, we present tests arguing that 
our measurement is free of several systematic effects, 
including effects arising from the data being obtained 
with different instruments; Ly/3, OVI, and associated 
absorption; and SLLS and DLAs. We also discuss how 
our continuum correction may itself suffer from some 
level of bias. However, our error bars thus far have 
not accounted for the systematic biases that may affect 
both the continuum and metal corrections. These biases 
are difficult to accurately quantify. A4oreover, they are 
most likely smoothly (thus affecting neighboring bins 
in with strong correlation) and monotonically varying 
with redshift. For the purpose of investigating whether 
the redshift evolution of T c ff shows evidence of a narrow 
feature (see, e.g., <J5|), these overall trends are mostly 
irrelevant and we will only consider the statistical errors 
estimated as above. Nonetheless, it is important to gain 
some idea of how uncertain the overall redshift evolution 
we report really is. We explore this question in this 
section. 

We choose to present t he calculation with the metal 
correction based on the ISchave et al.l (|2003f ) results as 
the fiducial model. The systematic error budget will 
account for the fact that the c orrect correction may 
be closer to that based on the iKirkman et al.1 (2005) 
results. 
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Fig. 8. — Comparison of our measurement of the total T e ff (red circles) with the inferred values including Lyct absorption only, r e £f iQ , 
with the metal contribution estimated based on the fits of Scliavc ct al. (2003) (blue squar es), who manua ll y ma sked metal lines in their 
spectra (left), and with the metal contribution statistically estimated based on the fits of Kirkman ct al. (2005) (green triangles). The 
metal-corrected data points have been slightly offset to the right for graphic clarity, but they are actually center ed on the same redshifts as 
the uncorrected r c ff measurement. All measurements have been corrected for continuum bias, as described in i\3A\ Numerical values are 
provided in Table l4l 



4.1. Continuum Bias Systematic Error 

In Appendix £|B.3[ we estimate the uncertainty on T e g 
arising from the continuum correction to |Ar c g|/r*™ c < 
0.0044(1 + z) . As this is a somewhat conservative 
upper bound estimate, we identify the la error to half 
this value: 

<Jr c!f , cent = 0.0022(1 + Z) 1 " 71 . (15) 

4.2. Metal Correction Systematic Error 

The systematic error on the metal correction is difficult 
to assess. Th e correction based o n the manual line iden- 
tifications bv lSchave et al.l |2003) is really a lower bound, 
as some lines are likely missed. This will necessarily oc- 
cur most frequently at high redshifts, where line blending 
is most important. On the other hand, metal lines co- 
incident with absorption already saturated by hydrogen 
have little to no imp act on the effective optical depth 
(e.g.. lKim et aLlfeOOTf ). so that it is unclear whether the 
larger fraction of missed lines at higher redshifts actually 
makes the estimated correction less acc urate. The sta- 
tistica l correction based on the results of Kirkman et alJ 
(2005) relies on a linear extrapolation from redward of 
the forest, which may not be exact. The two methods, 
however, are largely independent and a reasonably con- 
servative estimate of the error on the correction is ob- 
tained by taking the maximum deviation between the 
estimates from the two methods over the redshift inter- 
val probed. The deviation between the two corrections 
is defined as 

A^cff,a 1 7"cff,CK, Schaye ^eff, a, Kirkman] (16) 

Toff, a — T e ff T c ff iQ , Schaye — Teff 



where the subscripts 'Schaye' and 'Kirkman' refer to the 
measurements corrected for metal absorption as in equa- 
tions [13] and 1141 respectively, and T e ff is the continuum- 
corrected measurement prior to metal correction. It is 
largest in the lowest and highest redshifts bins, where it 
attains 60%. 9 We thus estimate the error contribution 
of the metal correction to the uncertainty on r e g as 

0V off , metal = 0.6(1 — T ff ,a,Schaye/T c ff) , (17) 

where the ratio T ff jQi s c h a ye/T e ff of the metal-corrected 
measurement to the value prior to correction is given by 
equation [TO] 

4.3. Total Error 

The total systematic error budget is obtained by also 
adding the systematic errors from the continuum and 
metal correction in quadrature, 

^Teff ,syst ^Yef^cont "T" ^*T e ff , metal ' (^^) 

The total uncertainty at each redshift is similarly ob- 
tained by adding the statistical contribution estimated 
as in i)3.3l and the systematic one in quadrature: 

^Voff.tot = a T att , stat + °V e « ,syst- (19) 

As discussed above, the systematic term is expected to 
be strongly correlated in neighboring redshift bins. The 
statistical error is estimated as in §3.3[ where it was 
simply denoted tr T( , ff . 

9 The one exception is the z = 2.1 bin for the Az = 0.1 mea- 
surement, for which the deviation is 73%. We regard this highest 
value as likely originating from a statistical fluctuation. 
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Fig. 9. — Illustration of the systematic error bu dget associated with the continuum and metal correct ions. The b lack point s sho w our 
r e ff measurements, corrected for continuum bias ( i|3,4H and with metals removed based on the results of Schayc ct al. (2003) f i|3,5l l, with 
the corresponding statistical errors calculated as in i|3.3l The dark gray bands include, in addition, the systematic contribution to the 
error budget from the uncertainty on the continuum correction. The light blue bands indicate the total error budget, with the uncertainty 
on the metal correction included as well. The bands serve to indicate the correlated nature of the systematic errors. The systematic 
contribution is at most comparable to the statistical error, and much smaller at both the low and high redshift ends. The left panel shows 
our measurement in Az = 0.2 bins and the right panel shows the measurement with Az = 0.1 bins. 

our continuum-corrected measurement in Az — 0.1 bins, 
with metals left in, is 

T 0ff = 0.0018(1 + zf m . (21) 

For the data with bins of width Az = 0.2 (Figure 
[T0|) . the z = 3.2 measured r e ff are seen to lie below the 
best-fit power laws at the 3a level, although the fits 
globally deviate from pure power laws only at < 2a. 
However, when the data are more finely binned with 
Az = 0.1 (Figure [TTj) . the reduced \ 2 for the power-law 
fits deteriorate substantially, and up to four consecutive 
points around z = 3.2 lie below the best-fit power laws. 
For example, our measurement corrected for continuum 
bias but including metals (first panel) has a probability 
of only 0.0007 of being described by the best-fit power 
law, with x 2 /dof = 48.0/(23 - 2) = 2.28. This is 
suggestive of a relatively narrow departure from a power 
law in the redshift evolution of r g near z — 3.2, which 
is partially smoothed out when the data are binned in 
redshift bins of width Az = 0.2. 

In order to investigate the possible departure from 
a power law around z — 3.2 in our mea surements, and 
to co mpare our results with those of iBernardi et al.l 
(2003), who found similar feature in their precision 
SDSS measurement, we also fit the same functional form 
consisting of a power law plus a Gaussian "bump" that 
these authors used, 



The results for both our measurements with Az = 0.1 
and Az = 0.2 bins are shown in Figure [5] The system- 
atic contribution is at most comparable to the statistical 
error, and much smaller at both the low and high 
redshift ends. 

5. FITS AND COMPARISONS WITH PREVIOUS WORK 
5.1. Fits 

In this section, we fit mathematical functions to our 
measurements of T c g versus redshift, and investigate 
how well each functional form fits the data. As we are 
interested in detecting potential narrow features in the 
redshift evolution of r e ff, we ignore the slowly varying 
systematic errors arising from the continuum and metal 
corrections, and instead consider only the statistical 
error bars estimated as in £13.31 We provide our results in 
graphical and tabular forms. Theoretical interpretation 
is left for separate work, with some possibilities outlined 
in S3 



The first functional form that we fit to our measure- 
ments is a power law, following previous studies that have 
found that the evolution of the Lya intergala ctic opacity 
was w ell described by such a function (e.g., iPress et al.l 
1993). Specifically, we solve for the parameters a and b 
which minimize the \ 2 between the model 



logr eff = a + 61og(l + z) 



(20) 



r cff = A(l + z) B + C cxp \ - 



and the measurements. These fits to our measurements, 
corrected for continuum bias and optionally for metal ab- 
sorption, are shown in the first rows of Figures [10] and 
ITTT for the data in redshift bins Az = 0.2 and Az = 0.1, 
respectively. The best-fit parameters and reduced \ 2 f° r 
each fit are tabulated in Table [5] In calculating x 2 /dof, 
we have subtracted the number of parameters fitted for, 
in this case two, from the number of data points to 
calculate the number of degrees of freedoms. For com- 
parison with previous results, the best-fit power law to 



[(1 + z)-D] 
2E 2 



(22) 



These fits are shown in the second rows of Figures [TU] 
and [IT] and the best-fit parameters and reduced \ 2 f° r 
each fit are also tabulated in Table [5] When calculating 
the reduced \ 2 f° r these fits, we have subtracted the 5 
parameters fitted for from the dof. For the data with 
bins of width Az = 0.2 (Figure fT0|) . the reduced % 2 are 
anomalously low (< 1), suggesting, as may have been 
expected, that the fitted model allows for too much 
flexibility given the number of data points. For the finer 
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Fig. 10. — Fits to our r e g measurements in redshift bins of width Az = 0.2. In all cases, the measurements have been corrected for 
continuum bias as in i]3.4l The top row shows power-law fits and the bottom row shows fits of a power law plus a Gaussian bump. In the 
first column, we fit directly to T ef j (including metal absorption). In the second and third column, we fit to T ef f a , our measurement with 
metals removed based on the fits of[Schayc et al. (20&3) and on the statistical technique of Ki rkman et al. (2005), respectively. For each 
fit, we give the reduced \ 2 value, where the number of parameters estimated from the data has been subtracted to obtain the effective 
number of degrees of freedom in each case. Best-fit parameters are given in Table Qj] 



Az = 0.1 bins (Figure [TTj) . the reduced. \ 2 are somewhat 
improved over the pure power-law fits. For example, 
it goes down from 2.28 to 1.81 (with corresponding 
p-values of 7 x 10 -4 and 0.012, respectively) for the 
measurement including metal absorption. However, even 
after the addition of the Gaussian bump near z — 3.2 to 
the power law, the data are not very well described by 
Equation (|2"2"|) . Indeed, the z = 2.3 data points exceed 
the best-fit models at the 3ct level and, although the 
evidence is weak owing to the large error bars, there is 
a hint of a roll over in the r e g evolution at z > 4. For 
instance, the measurement corrected for metal absorp- 
tion s tatistically based on the results of Kirkman et al. 
(120051) (last panel) yields x 2 /dof = 37.6/(23-5) = 2.01, 
corresponding to a p- value of 0.004. 

Inspection of Figure [TT] suggests that other func- 
tional forms may formally fit our data equally well 



as the power law+Gaussian model with a z « 3.2 
downward feature. This owes to the fact that the level 
of the error bars on our measurement allow for only a 
marginal detection of a departure from a pure power 
law on its basis alone. In particular, it does not clearly 
suggest the correct form of the deviation. The higher 
precision (but potentially affec ted by different sys - 
tematic effects) measurement of iBernardi et al.l (2003) 
did however clearly delineate a downward departure 
centered at z w 3.2, motivating us to fit the same 
functional form as these authors. Moreover, the redshift 
location and direction of the departure in this model 
are in qualitative agreement with the expectations of 
Hell reionization (see £JT|) . Nevertheless, we wish to 
emphasize that our data alone are unlikely to require 
this particular functional form. As there is no reason 
that T e ff should follow a power law even in the absence 
of a narrow feature, the parameters of the present fits 
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Fig. 11. — Fits to oui_T e ff measurements in redshift bins of width Az = 0.1. In all cases, the measurements have been corrected for 
continuum bias as in i]3.4l The top row shows power-law fits and the bottom row shows fits of a power law plus a Gaussian bump. In the 
first column, we fit directly to r e g (including metal absorption). In the second and third column, we fit to T ef f a , our measurement with 
metals removed based on the fits of[Schayc et al. (20Q3) and on the statistical technique of Ki rkman et al. (2005), respectively. For each 
fit, we give the reduced \ 2 value, where the number of parameters estimated from the data has been subtracted to obtain the effective 
number of degrees of freedom in each case. Best-fit parameters are given in Table [5] 

do not have particular physical significance. Rather, the red side of the Lya line, where there is no absorption 
our fits are merely a convenient device to test for the from the Lya forest. We begin by briefly describing the 
presence of a departure from a pure power law with an measurements of each class that overlap in redshift with 
implicit prior that the depart ure should be s imilar to ours, 
the one previously detected bv lBernardi et ail (|2003l ). 



5.2. Comparisons with Previous Work 

Further insight on the authenticity of the z ss 3.2 fea- 
ture suggested by the data is obtained by comparison 
with previous measurements. We compare with our mea- 
surement including absorption from metals, since many 
previous studies did not attempt to remove them. 

The existing T c g measurements fall in two broad 
classes: (1) measurements from high-resolution spectra, 
for which the continuum level of each spectrum is 
directly estimated on an individual basis, and (2) 
measurements from low-resolution spectra, with the 
continuum level estimated from an extrapolation from 



5.2.1. Descriptions of Existing r e a Measurements 

For the direct measurements from high-resolution spec- 
tra, we focus on those that were obtained from > 20 
spectra and for which the evolution with redshift of T e s 
h as been estim a ted. 

iSchaye et al.l (|2003h manually estimated the contin- 
uum level of 14 UVES and 5 HIRES spectra, in a 
manner similar to ours, with quasar emission redshifts 
2.145 < z em < 4.558. As discussed in §3.51 they re- 
peated their measurement with and without removing 
metal lines; here we consider their measurement with 
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TABLE 5 

Parameters for the best-fit power-law and power-law+Gaussian models to our measurements of r e g, 

WITH AND WITHOUT METAL CORRECTIONS, AND IN REDSHIFT BINS OF WIDTH Az = 0.2 AND Az = 0.1 



Power law a Power law+Gaussian b 

b x 2 /dof c A B C D E x 2 /dof d 

Az = 0.2 



Including metals -2.706 3.880 1.74 0.00196 3.892 -0.0708 4.207 0.101 0.54 

Schaye metal correction -2.853 4.057 1.48 0.00141 4.062 -0.0651 4.207 0.101 0.47 
Kirkman metal correction -2.910 4.164 1.63 0.00115 4.222 -0.0673 4.218 0.0900 0.86 

Az = 0.1 

Including metals -2.734 5T924 2~28 0.00153 4U60 -0.0969 OB7 0.0769 1781 

Schaye metal correction -2.876 4.094 1.91 0.00111 4.225 -0.0891 4.267 0.0766 1.49 
Kirkman metal correction -2.927 4.192 2.31 0.000876 4.403 -0.0930 4.269 0.0726 2.09 
a Model defined in Equation Il20jl . b Model defined in Equation J22jl. c We have subtracted the two parameters 
fitted for from the number of degrees of freedom used to calculated the reduced x 2 - d We have subtracted the five 
parameters fitted for from the number of degrees of freedom used to calculated the reduced x 2 ■ 



metal absorption included. They did not apply an ex- 
p licit corre c tion fo r redshift-dependent continuum bias. 

iSongailal (|2004f ) studied 25 HIRES spectra with emis- 
sion redshifts 2.31 < z em < 4.10 and 25 ESI spectra 
with emission redshifts 4.17 < z em < 6.39. She manu- 
ally estimated the continuum levels for the quasars with 
Zem < 4.5, but extrapolated from the red side of the 
Lya line for the higher-redshift ones. For the highest- 
redshift quasars in her sample, she used a fixed power-law 
spectral index a = —1.25 (/„ oc v~ a ) and set the nor- 
malization using the region of the spectra near the rest- 
frame wavelength 1350 A. She did not apply a redshift- 
d ependent contin u um co rrection. 

Kirkma n et al.l (|2005f ) measured the flux decrement 
DA at 1.6 < z Lya < 3.2 using 24 HIRES spectra. 

Following the earlier work of iTvtler et al.1 (|2004al ) at 
z w 1.9, they had a team of undergraduate students con- 
tinuum fit both their actual spectra and mock spectra 
with matching properties. They estimated a small bias 
as a function of DA, which they used to correct their 
measurement. The maximum correction was w 1%. At 
2 = 3. our continuum bias correction w 4% is larger than 
their maximum correction, althoug h our results a re not 
necessarily inconsistent. In fact, iKirkman et alj (2005) 
did not estimate a correction as a function of redshift 
and most of their data is at z < 3. Instead of using 
a co smological simu l ation to generate their mock spec- 
tra, IKirkman et al.l (|2005h randomly placed Voigt pro- 
files, f ollowing the line distributions of lKirkman fc Tytlerl 
(1997) and with random redshifts. These mock spec- 
tra thus neglect correlations between the absorption 
lines and do not rep resent the true cosmic structure. 
IKirkman et al.l (|2005l ) find that these mock spectra con- 
tain w 5% more total absorption than their real ones. 
These differences may thus explain discrepancies that 
may (or may not) exist between their continuum cor- 
rectio n and ours. As discussed in H3.51 Kirkma n et al.l 
( 2005) statistically estimated and subtracted the absorp- 
tion arising f rom metals in t he Ly a forest using line lists 
p ublished bylSargent et all (|1988h . 

iBecker et al.l (|2007|) most recently measured the prob- 
ability density function (PDF) of the transmitted flux F 
in the Lya forest of 55 HIRES spectra, with absorption 
redshifts spanning 1.7 < z^ a < 5.8. At z^ a < 5.4, 

they used a manual spline fit to estimate the quasar 
continua. They noted that the transmitted regions at 



z Lya — ^ rarely, if ever, reached the continuum, and 
so their fits were of very low order for these redshifts. 
They used a power-law fit to the continuum on the red 
side of the Lya line with a = 0.5 for quasars with 
Zem > 5.7. Beck er~et al.l ()2007f ) removed all the metal 
lines that could be identified either by damped Lya ab- 
sorption or from multiple metal lines at the same redshift. 
They did not, however, mask weak lines found in the for- 
est without counterparts redward of Lya emission, with 
the motivation that doing so would preferentially discard 
pixels with low Lya optical depth (where the metal lines 
can be seen), introducing a potentially larger bias in the 
PDF than the one incurred by leaving the contaminated 
pixels in the sample . From lognormal fi ts to the observed 
optical depth PDF. IBecker et all (|2007f ) inferred the cor- 
responding r ff redshift evolution, but did not attempt to 
estimate the uncertainty in this derived quantity. They 
did not apply any continuum correction in this calcu- 
lation, but note d that if the functiona l form of the gas 
density PDF of iMiralda-Escude et all (|2000f ) is instead 
assumed, then a continuum correction similar in magni- 
tude to ours must be applied in order to obtain a good 
fit. 

Other works that have directly estimated the transmis- 
sion of the Lya forest from hi gh-resolution spectra, bu t 
using sma l ler da t a sets, include I Steidel &: SargentJ (I I 987), 
IHu et al.l (|1995f). iRauch et all (|1997l ). iMcDonald et all 
(|2000f ). and lKim et all (|2001l ). To compare our study to 
previous analyses of the Lya forest from high-resolution 
spectra, we note that the two previous studies which 
included the largest samples of q uasars over the red- 
shift range we cove r are those of ISongailal (|2004l ) and 
IBecker et al.l (l2007f) . Thirty of the quasars analyzed 
by Songailar(|2004f ) have a median Lya absorption red- 
shift between 2 and 4, where our measurement is most 
sensitive (Figure [H). Thirty - nine of the quasars ana- 
lyzed by IBecker et al.1 (|2007l) contribute a point with 
mean Lya absorption redshift in this interval. Our mea- 
surement from 86 quasars thus represents an improve- 
ment over previous analyses of the Lya forest from high- 
resolution spectra at 2 < z < 4 of a factor of about two 
in the size of the data set alone. 

The measurements using low-resolution spectra, with 
extrapolation of the continuum fr om the red side o f 
the Lya line, were pioneered by iPress et al.l (l993). 
These authors analyze d 29 quasars from the sample of 
ISchneider et alj (| 19911 ) with a resolution of 25 A. They 
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extrapolated a continuum of the form C1/2A 1 / 2 + C\X, 
where C1/2 and C\ are parameters that are fit for, ac- 
counting for the presence of emission lines redward of 
Lya emission. They ignored possible features, such as 
emission lines, in the Lya forest region. They fitted a 
power law of the form r c ff = A(l + z) 7+1 to their data 
over the redshift range 2.5 < < 4.3, obtaining 

A = 0.0175 - 0.0 567 ± 0.0002 and 7 = 2.46 ± 0.37. 

iBernardi et a l. (2003) analyzed 1061 spectra from the 
SDSS and from those measured the evolution of T e g over 
the redshift range 2.5 < z^ a < 4. They fit a continuum 

consisting of a power law plus three emission lines (Lya, 
1073 A, and 1123 A) in the Lya forest region, fiducially 
normalizing each spectrum by its flux in the rest-frame 
wavelength range 1450-1460 A. They found that the T e ff 
evolution was well-fitted by a power law with 7 = 3.8 ± 
0.2, except for a clear feature taking the form of a "dip" 

n ear z = 3.2. 

iMcDonald et al. (2006) performed a similar analysis on 
an expanded sample of 3035 SDSS spectra, but using a 
principal components analysis to model the quasar con- 
tinuum shape, and normalizing their measurement us- 
ing the rest-frame wavelength range 1268-1380 A. Their 
measuremen t of the mean transmissi on (F) (z) is given in 
Figure 20 of IMcDonald et all (|2005l ) and does not show 
evidence for a pronounced feature near z — 3.2, although 
these authors warn that their data points are correlated, 
that their error bars may not be accurate, and therefore 
that their measurement should be used for "qualitative 
use only" . 

The contribution of metals to absorption in the 
Lya forest was not explicitly sub tr acted in the 

iPress et al.1 dl993j ^ IBernardi et~ai1 (|2003l ). and 

McDon ald et alJ ( 2005T ) measurements. However, 



extrapolating the continuum from redward of Lya emis- 
sion presumably accounts for absorption by metallic 
transitions with rest-frame wavelength longer than 
Lya to some level. 



5.2.2. Comparisons with Other Direct Measurements 

The published z > 3 points for the direct measure- 
ments, shown in the left panel of Figure [121 arc in better 
agreement with our measurement prior to correc t ion fo r 
continuum bias , as e xpected since iSchaye et al.1 l(200l . 
ISongailal (|2004l ). and IBecker et all (120071 ) did not apply 



similar corrections. Note that none of the published data 
points for the direct measurements exceeds our best-fit 
power law near z — 3.2 (even prior to continuum correc- 
tion), consistent with the apparent feature in at this red- 
shift in our meas urement. We note, i n passing, that the 
error bars on the lSchave et alJ (|2003l ) points appear un- 
derestimated, based on comparing them with the scatter 
between th e individual da ta points. Similarly, the error 
bars on the ISongailal (|2004l ) points appear overes timated. 

The r e ff curve inferred bv lBecker et al.l (|2007h appears 
to qualitatively differ from our measurement. These 
authors did not, however, attempt to directly measure 
(F) or r c ff. Instead, they assumed that the redshift 
dependence of the logarithmic mean /i and standard 
deviation a of their lognormal fits to the flux PDF 
was well approximated by a linear function in 1 + z. 



From these best- fit linear functions, they calculated 
(F)(z), and hence r e g(z). The shape of their inferred 
r e g curve reflects the assumption of linearity of /i 
and a in 1 + z, which may not be exact, especially if 
the evolution of r e s has a feature near z = 3. Their 
measurement is also expected to lie somewhat below 
ours, since these authors did not apply a continuum 
correction t o the measur e ment used to produce their 
r ff curve. IBecker et al.1 ((2007) did not quantify the 
uncertainty on this derived quantity, but it is expected 
to be significantly larger than the errors on our mea- 
surement (G. Becker, private communication), in part 
because of the smaller number of spectra contributing 
to the redshift range 2 < z^y a < 4.2 in their data 
set. A small fraction of the difference be tween our 
measurement and the IBecker et al.1 ((2007) curve is 
attributable to the fact that these authors removed 
most of the metal lines in their measurement, whereas 
we left them in for the measurement plotted in Figure [T2l 



5.2.3. Comparisons with Continuum Extrapolation 
Measurements 

Our measurement corrected for continuum bias is 
seen to trace the S DSS statistical measurement of 
IBernardi et al.1 (|2003f ). who first claimed detection of a 
feature at z = 3.2, remarkably well. In Figure Q21 we 
quantitatively compare our meas urement of r e ff correcte d 
for continuum bias wit h that of | Bernardi et al.l ([20031 ) . 
To do so, we binned the IBernardi et aTl J2003) points in 
the same Az = 0.1 redshift bins as our measurement. 
The error bar on the binned points is estimated as the 
quadratic sum of the data points that fall in a partic- 
ular bin, divided by the square root of the number of 
points in the bin. These error bars are slight under- 
estimates, sinc e nei ghboring data points published by 
IBernardi et al.l (|2003l ) are slightly correlated. The true 
error bars would only render our conclusion stronger. We 
have a dded a constant vertica l offset in logarithmic space 
to the IBernardi et alJ (|2003T ) points that minimizes the 
X 2 between the two data sets (AlogT c g = 0.011), since 
their statistical measurement may not be accurately nor- 
malized 10 . Accordingly, we have subtracted one from the 
number of degrees of freedom used to calculate the re- 
duced x 2 - The resulting % 2 /rfo/ = 16.4/(14 - 1) = 1.26 
indicates that the two data sets agree very well, with 
a p-value of 23%. T he z « 3.2 feature discovered by 
IBernardi et al.l (|2003f ) has sometimes been suspected of 
being an artifact of the SDSS spectrograph, which has 
a dichroic optical element splitting the beam into blue 
and red arms near the observed wavelength correspond- 
ing to Lya at z ~ 3. The ESI and HIRES spectrographs 
on Keck, with which 60 out of the 86 spectra we ana- 
lyzed were obtained, consist of a single arm and there- 
fore do not suffer from the same potential systematic ef- 
fect. Moreover, in this work, we locally estimate the con- 
tinuum level of each spectrum and so our measurement 



10 T he normalization on the r c g measurement of Bcrnardi et al. 
(2003) was obtained by normalizing each spectrum by its flux in 
the rest-wavelength range 1450-1470 A. As these authors show, 
their absolute r e fj measurement depends somewhat on the rest- 
wavelength range used for the normalization. 
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Fig. 12. — Comparisons of our r c ff measurements with published results. The left panel compares measurements in which the quasar 
continuum levels were estimated directly from high-resolution, high-signal-to-noise spectra, as in this work. The right panel shows estimates 
based extrapolating the quasar continuum from redward of Lya emi ssion in low-resolution spectra. The solid black curves show the best- fit 
power law to our measurement corrected for continuum bias ( i|3.4l l. but including metal absorption. Residuals are given with respect to 
these curves. The dashed black curves indicate the best-fit power law to our measurement prior to correction for continu um bias. The 
publis h ed z > 3 points are in better agreement wi th our measurement prior to correction for continuum bias, as expected since ISchave et ahl 
(2003), Songaila (2004), and Becker et al. (2007) did not apply similar corrections. Note that none of the published data points for the 
direct measurements exceeds our best-fit power law near z = 3.2 (even prior to continuum correction), consistent with the apparent feature 
at this redshift in our measurement. Our measurement corrected for continuum bias is seen to trace the SDSS statistical measurement of 
[Bernardi et a l. (2003), who have first clai med detection of a feat ure at z = 3.2, remarkably well. We make this comparison quantitative 
in Figure [13] The SDSS measurement of McDonald et al. (2005) deviates significantly from ours at z > 3, but these authors warn that 
their points are co rrelated and that their error bars may be underestimated. Our measurement agrees well with the best-fit power law of 
Press et al. (1993), although the uncertainty on the latter is large. 



should be largely unaffected by flux calibration problems, 
since it is only a function of t he normalized transmiss ion. 

The SDSS measurement of lMcDonald et all (|2005f ) de- 
viates significantly from ours from below at z > 3, al- 
though these authors warn that their points are corre- 
lated and that thei r error bars may be un derestimated. 
It is of note that the lMcDonald et al.l (|2005f ) r e ff measure- 
ment lies below even our raw measurement uncorrected 
for continuum bias, which should be a rather firm lower 
bound to what r e ff can be. Part of the discrepancy may 
be attributed to the fact that absorption arising from 
metals present redward of th e Lya line is impl i citly s ub- 
tracted to some extent in the lMcDonald et all ()2005l ) es- 
timate. However, the discrepancy in the highest redshift 
z = 4.2 bin is ~ 30%, while our estimate of the metal 
absorption at this redshift is only 2 — 5% (Sj33|. Recon- 
ciling the measurements thus appears to require errors on 
the continuum and metal corrections much larger than 
the systematic uncertainties estimated in <21 indicating 
a genuine discrepancy 

The best-fit power law of iPress et al.1 £[993) agrees 
with ours within the uncertainties, although the error 
on the former is large. 



6. DISCUSSION 



We have measured the evolution of the Lya effec- 
tive optical depth over the redshift range 2 < z < 4.2 
from a sample of 86 high-resolution, high-signal-to-noise 
quasar spectra obtained with Keck/ESI, Keck/HIRES, 
and Magellan/MIKE, paying particular attention to ro- 
bust error estimation and potential systematic effects. 
This represents an improvement over previous analyses 
of the Lya forest from high- resolution spectra in this red- 
shift interval of a factor of two in the size of the data set 
alone. 

Using mock spectra, we have calculated the redshift- 
dependent bias on our measurement arising from the in- 
creasing cosmological density - and hence absorption - 
with redshift ( §3.4j) . This correction is important at red- 
shifts z > 3, with a magnitude w 12% at z = 4. Previous 
direct r c ff measurements, in which the quasar continua 
were estimated by fitting the transmission peaks in the 
Lya forest, at these redshifts generally did not apply 
such a correction and are thus likely to be biased low. 
This effect docs not affect measurements based on ex- 
trapolating the continuum from redward of the Lya line, 
as usually done in measurements based on lower res- 
olution data for which direct local continuum estima- 
tion i s impossible (e.g.. IPress et al.lll993l : iBernardi et al.1 
120031: iMcDonald et al.ll2005fh Such extrapolations, how- 
ever, likely miss local variations in the continuum from 
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Fig. 13. — Comparison of our m easurement of T e g corr ected for continuum bias (see 33.41 red circles) and of the Bernardi et al. (2003) 
measurement (blue squares). The IBernardi et al. (2003) points were binned in the same Az = 0.1 redshift bins as our measurement. 
The blue data points are cen tered on the same red shifts as the red ones, but have been slightly offset to the right for graphical clarity. 
The error bars on the binned Bernardi ct al. (2003) points arc slightly underestimated, since their neighboring published data points are 
slightly correlated. We have added a constant vertical offset in logarithmic space to the Bernardi et al. (2003) points that minimizes the 
X 2 between the two data sets (AlogT c g = 0.011), since their statistical measurement may not be accurately normalized. Accordingly, we 
have subtracted one from the number of degrees of freedom used to calculate the reduced x 2 - The resulting x 2 /^o/ = 16.4/ (14 — 1) = 1.26 
indicates that the two data sets agree very well, with a p-value of 23%. 



quasar to quasar, making local estimation advantageous 
for high-resolution data. As discussed in Appendix [El 
the continuum correction estimated in this work is itself 
subject at some level to systematic uncertainty, arising 
principally from uncertainties in the cosmological hydro- 
gen photoionization rate, T bkg . This systematic uncer- 
tainty could in principle be eliminated by iteratively re- 
estimating the continuum correction and T bk9 from the 
data itself until convergence is attained. We defer this 
involved procedure to future work. 

We have also estimated the level of absorption in the 
Lya forest arising from metals, based on results in the 
literature on both direct and statistical metal removal, 
finding that this contribution is rs 6 — 9% at z = 3 and 
decreases monotonically with redshift f £)3.5p . 

The high precision of our measurement (attaining 3% 
in redshift bins with width Az = 0.2 around z = 3), 
indicates significant departures from the best-fit power 
law evolution (r c ff = 0.0018(1 + z) 3 92 for our measure- 
ment corrected for continuum bias, but including metals; 
Table [5]), particularly near z = 3.2, where there is evi- 
dence of a dip (SJSJ). Our measurement, with the feature 
near z = 3.2, is in excellent agreement with the preci- 
sion m easurement from SDSS spectra of IBernardi et all 
(|2003l) (^5.21). who used an independent approach. 

iTheuns et al.l (l2002aTl interpreted the feature seen by 
IBernardi et al.l (|2003f ) as a signature of Hell reionization. 
Although the feature is indeed in qualitative agreement 
with what might be expected from a temperature jump 
in the IGM owing to Hell reionization, degeneracies ex- 



ist. Neglecting redshift-space distortions, the Lya optical 
depth at any point in the forest is given by 



ne2 fLya 1 R(T)n ml n e 



Ybkg 



(23) 



where we have m ade the approximation of photoioniza- 
tion equilibrium (|Gunn fc Petersonl 1 19651 ). Here, e is 
the electron charge, m e is the electron mass, /t 



yo 



is 



the Lya oscillator strength, ^Lyoj ^ s ^ ne ^y a frequency, 
R(T) oc T-°- 7 is the hydrogen recombination coefficient 
(e.g.. lOsterbrock &; Ferlandll2006f ) . nmi is the ionized hy- 
drogen number density, n e is the free electron number 
density, and T bk9 is the hydrogen photoionization rate. 
From Equation (f2"3"]l . we see that at least three things 
could produce a feature in the redshift evolution of the 
optical depth: 

1. A change in temperature: an upward jump in T 
will result in a downward jump in t; 

2. a change in the free electron number density: an 
upward jump in n e will result in an upward jump 
in r; 

3. a change in the ionizing background: an upward 
change in T bkg will result in a downward change in 

T. 

Both an upward temperature change and an increased 
free electron number density could arise from the reion- 
ization of Hell, although the latter effect is limited to 
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7% if Hel has already been reionized (e.g., by stars si- 
multaneous to HI reionization), for a cosmic He mass 
fraction Y = 0.25. A temperature jump of a fac- 
tor of two (jRicotti et al.ll2000t ISchave et al.ll2003l ) could 
suppress r by up to 2" 7 « 40%, although redshift- 
space distortions coul d mitigate this effect at some level 
(jTheuns et all l2002a| ) . However, even if Hell reioniza- 
tion is not happening and the IGM temperature is con- 
stant, r e fi could well display a feature if the ionizing 
background is evolving suitably, for example if it is tran- 
sitioning from being dominated by stars to being domi- 
nated by quasars. F or instance, empirical studies (e.g., 
Hopk ins et al.l l2007f ) indicate that the contribution to 
the ionizing background from quasars peaked at Z f» 2, 
with small uncertainty, but simple theoretical arguments 
(jHernquist fc S pringcl 2003) imply that the contribution 
from stars likely peaked much earlier. We must note 
that the effective optical depth we have measured is not 
directly related to the local optical depth as given by 
equation[23j but in fact depends on the PDF of the trans- 
mitted flux through its first moment (eq. 

In addition to the degeneracies outlined above, the the- 
oretical interpretation of the feature we find evidence for 
near z = 3.2 poses significant challenges. Why does the 
feature appear so narrow? If it is due to a tempera- 
ture jump of a factor of two from Hell reionization and 
the IGM cools mainly adiabatically, then the Universe 
must expand by a factor of \[2 in order for the IGM 
temperature to go back to its pre-jump value. Yet, the 
feature appears to have a much smaller width (Az ~ 0.5 
at z ~ 3.2; e.g., Fig. [IT]). Cosmic variance may also 
be expected to smooth the observed feature, which is 
an average over many sightlines, although it is not clear 
whether this would act to increase or decrease its width. 
Moreover, our estimate of the duration of the event may 
be misled by our perception of r e g evolution as a power 
law plus a narrow downward feature. Indeed, even in the 
absence of an event such as Hell reionization near z = 3, 
there is apparently no good reason to expect r c g to evolve 
as a pure power law, so that the baseline against which 
we are comparing the feature may not be correct. A fea- 
ture in the redshift evolution of the mean transmission 
of the Lya forest may also be expected to impact its flux 
power spectrum. Precision meas urements of the large - 
scale flux power from SDSS (e.g.. lMcDonald et al] l2006) 
apparently do not show a corresponding feature. Prelim- 
inary calculations, however, suggest that the large-scale 
power need not be affected by a feature in the mean 
transmission of the forest, if it arises from a change in 
the temperature of the gas, since a change in the tem- 
perature modifies the range of gas densities to which the 
power spectrum is sensitive in a way that tends to coun- 
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teract the effect of the change in mean transmission. 

We regard the questions of whether Hell reionization 
is occurring near z — 3.2, whether the detection of a dip 
in the evolution of the Lya effective optical depth near 
this redshift could be a signature of this process, or what 
are the other implications of our precision measurement 
of the evolution of Lya effective optical depth, as rich 
and complex. We will investigate these questions theo- 
retically in separate work. 

Observationally, additional clues on the evolution of 
the physical state of the IGM at 2 < z < 4 could be pro- 
vided in the immediate or near future by measurements 
of higher-order (/3, 7, ...) Lyman transitions, which probe 
different density regimes. Precise measurements of the 
small-scale flux power spectrum of the forest, which is 
perhaps most sensitive to the temperature of the IGM, 
with fine redshift sampling may offer our best oppor- 
tunity to break the degeneracy between a temperature 
change and a pure evolution of the hydrogen photoioniz- 
ing background. 
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APPENDIX 
TESTS FOR SYSTEMATIC EFFECTS 

We present in this section tests that we performed to ensure that our measurement is free of significant systematic 
errors. 

Different Instruments 

As detailed in S}21 the spectra we analyzed in this work were obtained with the HIRES and ESI spectrographs on 
Keck, and with the MIKE spectrograph on Magellan. To verify that the data obtained with the different instruments 
give consistent results, we have repeated our raw r c g measurement (uncorrected for continuum bias and with metals 
left in), separately using only data from Keck (HIRES+ESI), then only data from Magellan (MIKE). We have chosen 
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Fig. 14. — Comparison of T c ff (uncorrected for continuum bias and with metals left in) as measured from the Keck (HIRES+ESI; red 
circles) and Magellan (MIKE; blue squares) spectra only. The blue data points are centered on the same redshifts as the red ones, but 
have been slightly offset to the right for graphical clarity. The reduced \ 2 = 0.90 between the two measurements indicate that the data 
from the different observatories give perfectly consistent results. 



to group the HIRES and ESI data into Keck data because the small number of HIRES spectra (16, distributed as 
illustrated in Figure [2]) makes it difficult to calculate reliable error bars on the measurement from these spectra alone. 
The result of this comparison is shown in Figure [Ml The reduced \ 2 = 0.90 between the two measurements indicate 
that the data from the different observatories give perfectly consistent results. Figure[S]in i)3.4| provides a similar test, 
combining the data from different redshift bins by normalizing them according to the mean transmission in the bins. 

Ly(3-OVI Emission and Associated Absorption 

Some authors (e.g. JBecker et al.ll2007fl have previously masked a region redward of Ly/3 in their Lya forest analyses 
to avoid potential problems arising from the presence of the Ly/3 emission line itself, as well as from OVI emission 
at 1038 A. In this work, it was unclear that these emission lines were problematic. We thus chose to keep all 
Lya forest wavelengths redward of Ly/3, except for the proximity regions of the quasars, in our analysis ( EI3.2.1j) . We 
repeated our raw measurement (with metals left and uncorrected for continuum bias) masking rest-frame wavelengths 
< 1050 A and found no significant difference, as shown in the top left panel of Figure [T"5l 



Quasar spectra occasionally show evidence for metal absorption associated with the host galaxies. We tested whether 
this could affect our measurement by masking metal lines (as for those associated with SLLS and DLAs, as described 
in i j3.2|) with redshift within 3000 km s _1 of our quasars. The width of these masks was chosen to account for the large 
typical uncertainties in quasar redshift determinations (for related discussions on quasar redshift uncertainties, see e.g. , 
GaskeUll983lTvtler fc FaiJ|1992t IVanden Berk et al.ll200lt iRichards et"alll2002t iHennawi et"aT]|2006t IShen et al.ll2007l ; 
Faucher-Giguere et al][2007l ). As the top right panel of Figure ITSI shows, masking this potential associated absorption 
also has no significant effect on our measurement. 

Super Lyman- Limit Systems and Damped Lya Absorbers 

Many of the lines of sight analyzed in this work were originally selected to contain super Lyman-limit systems 
and damped Lya absorbers. In order to make an unbiased measurement of r e s, we have masked the systems with 
log ./Vhi > 19.0, which individually contribute relatively large equivalent widths to the absorption, as described in £)3.2l 
In this section, we justify our claim that the resulting measurement should be nearly unbiased and test for systematic 
effects that may be introduced by the masks. 

We first show that in a random sample of sightlines, the contribution of systems with log /Vhi > 19.0 to r c g is small. 
Thus, if these systems are masked from a biased sample, the unmasked regions should have a r e ff close to that of a 
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Fig. 15. — Tests for systematic effects associated with choices of analysis parameters. In each panel, the red circles show our raw r e ff 
measurement (with metals left in and uncorrected for continuum bias) and the blue squares show the corresponding measurement for a 
variation of the analysis parameters. Top left (Ly/3 and OVI emission): Rest-frame wavelengths < 1050 A masked. Top right (associated 
absorption): Metal lines with redshift within 3000 km s — 1 of the quasars masked. Bottom left (width of SLLS and DLA masks): Masks 
around SLLS and DLAs twice as large as the fiducial width (SLLS or DLA transmission < 99%). Bottom right (masks on metal lines 
associated with SLLS and DLAs): Not masking metal lines associated with SLLS and DLAs. The blue points are centered on the same 
redshifts as the red ones, but have been slightly offset to the right for graphical clarity. In each panel, the two measurements agree within 
la over the entire redshift range covered, indicating that our results are robust to the choices of analysis parameters. 



random sample. 11 

We can estimate the contribution of SLLS and DLAs to the opacity of the Lyct forest from the observed column 
density distribution. Specifically, the average effective optical depth owing to systems with column density in the range 
[N min , N max ] can be written as 

(r eff ) = t / dN m / dbf(N m , b)W(N m , 6), (Al) 

A Lya JN min J 

where /(Am, b) is the joint PDF of the column density and the Doppler ^-parameter, andW(Njn, b) is the equivalent 
width of an absorber with column density Am and a given 6-parameter (e.g.. lZuolll993l ). We set logA mm = 19.0 
and log N max = 22.0. This maximum was chosen so that systems of column density this high are too rare to make a 
differ ence even in our biased sample, owing to the cutoff that is observed in the distribution of DLA column densities 
(e.g.. lSchavdl2lfflll ). We assume for simplicity that Am and b are statistically independent, and that the ^-parameter 
distribution is a S- function centered on b = 30 km s . As expected, since the equivalent width depends only weakly 
on the 6-parameter for the large column densities of interest, setting b = 20 km s -1 or b = 40 k m s _1 instead makes 
no practical difference. For the column density distribution, we use the H aardt fc Madaul (|1996f ) parameterization, in 
which /(A H i) oc A HI 15 . 

Figure [TBI shows the contribution of discrete absorbers with column density log Am > 19.0 to the total effective 
optical depth of the Lya forest as a function of redshift in an unbiased sample of sightlines. This contribution is 
fractionally small (< 2%) over the redshift range covered by our t c q measurement. Therefore, if these systems are 
completely masked in a biased sample, r c g is underestimated by this small amount with respect to an unbiased sample. 



11 There is an implicit assumption in this argument that the masked regions are uncorrelated with the unmasked ones. This assumption 
is satisfied, since the SLLS a nd D LA masks typically have widths of tens of proper Mpc, compared to a correlation length of about 1 proper 
Mpc in the Lya forest (see ^13.31 for a discussion of this correlation length). 
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Fig. 16. — Contribution of discrete absorbers with column density logA^i > 19.0 to the total effective optical depth of the Lya forest 
as a function of redshift in an unbiased sample of sightlines. This contribution is fractionally small (< 2%) over the redshift range covered 
by our r c fj measurement. 

Although the above argument is in principle sufficient, we have also verified that our results are insensitive to the 
width of our SLLS and DLA masks, thus confirming that no artificial trends are introduced by them. Specifically, we 
compare in the bottom left panel of Figure [15] our uncorrected measurement with the same measurement but with 
masks twice as large, and find the two measurements agree within la over the entire redshift range covered. The 
bottom right panel of Figure [15] shows that the measurement is also insensitive to masking the metal lines associated 
with the SLLS and the DLAs. 

CONTINUUM ERROR DETAILS 

In this section, we estimate the errors made when continuum fitting quasar spectra manually, as done to measure 
the evolution of the Lya effective optical depth in this paper. 

Variance Between Continuum Fitters 

First, we ask how much do continuum fits of the same spectra vary when they are made by different individuals? 
A sub-sample of 31 of our ESI spectra were independently continuum fitted by two of us, J. X. Prochaska and C.-A. 
Faucher-Giguere. The continuum fits were found to agree surprisingly well between the two continuum fitters, as 
illustrated in Figure [TT] The r e ff values from the different continuum fits agree within 2% on average, and everywhere 
within 4%, for redshift bins of width Az — 0.2. The discrepancies are generally much smaller than the la errors on 
the measurements, indicating that variance in the continuum fits is a negligible source of error. 



Systematic Bias 

Having determined that the scatter between different continuum fits of the same spectra is small, we turn to the 
question of whether these may be systematically biased with respect to the true quasar continuum. Of course, the true 
continuum level is unknown in actual quasar spectra. We therefore make use of mock spectra, for which the continuum 
level is known by construction. 

We generate transmission fields F(X) from a cosmological dark matter simulation, as described in detail in 
iFaucher-Giguere et al.l (|2007| ). Briefly, the baryons are ass umed to trace the dar k matter, except for Jeans smoothing, 
and obey a simple power-law temperature-density relation (|Hui fc Gnedin|[l997f ). The redshift-space distortions owing 
to thermal broadening and peculiar velocities are explicitly modeled. The simulation box we use has 1024 3 particles with 
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Fig. 17. — Comparison of r e g as inferred from a subset of 31 ESI spectra, where the spectra have been continuum fitted independently 
by CAFG (red circles) and by JXP (blue squares). The blue data points are centered on the same redshifts as the red ones, but have 
been slightly offset to the right for graphical clarity. The r e g values from the different continuum fits agree within 2% on average, and 
everywhere within 4%. The discrepancies are generally much smaller than the la errors on the measurements, indicating that variance in 
the continuum fits is a negligible source of error. 

side length 130 h" 1 co-moving Mpc. Fo r the hydrogen photoionizing background, T bkg , we li nearly interpo l ate b etween 
the values given bv lBolton et al.1 (j2005l) at 2 < z < 4; at z = 4.5, we interpolate between the lBolton et al.l (j2005) z = 4 
value and the z = 5 value given in lMeiksin fc White! (|2004l) . This results in T bk 3 = (1.3, 1.1, 0.9, 0.95, 1.0, 0.655) x 10~ 12 
s- 1 at z= (2,2.5,3,3.5,4,4.5). 

As real quasar spectra are not featureless, but display a number of weak emission lines in Lya forest, we multiply our 
normalized transmissio n fields by a composite quasar spectrum made from a sample of 657 quasars from the FIRST 
Bright Quasar Survey (|Brotherton et al.ll2001h . The mock spectra are degraded to realistic resolution, pixel size, and 
signal-to-noise. 

Careful manual continuum fitting is a very time-consuming procedure, and so it is impossible to explore the full 
redshift, resolution, and signal-to-noise parameter space covered by our actual spectra. The tests of ^3.41 support 
that there is little, if any, trend of continuum error with noise and resolution for our high-quality data. We thus 
concentrate on the particular case of ESI spectra, for which we have the most data at redshifts z > 3, where we 
expect the continuum error to be largest, owing to the increasing cosmological matter density. The ESI spectra are 
also those with the coarsest resolution used in this work. Specifically, we set FWHM = 40 km s -1 , 11 km s _1 pixels, 
and S/N = 20 per pixel. We generated 10 different mock spectra with these properties at redshift intervals Az = 0.5 
from z — 2 to z = 4.5, and continuum fitted each of them with the same procedure with which we fitted the continua 
of the real spectra. As a check, we also continuum fitted a few mock HIRES spectra (FWHM = 6 km s" 1 , 2 km s _1 
pixels, and S/N — 30 per pixel), corresponding to the highest resolution spectra in our sample, and confirmed that the 
continuum correction (see below) for these was approximately the same as for the ESI spectra at z = 4. The true mock 
continuum, while known by construction, was hidden during the continuum fitting. In generating the mock spectra, 
we wrapped around our simulation box 6 times for each spectrum, for a total path length of 6 x 130 hr 1 comoving 
Mpc. Each mock spectrum thus contains more Lya forest data than an actual spectrum would between 1025 A and 
1216 A, by a factor w 2 at z — 3. We intentionally did not model the cosmological evolution of the Lya forest in 
a given quasar spectrum versus observed wavelength. For example, the entire Lya forest of one of our z = 4 mock 
spectra arises from a z — 4 density field subject to a z = 4 photoionizing background. Therefore, the mean continuum 
errors that we calculate at a given redshift are valid for this exact redshift. 

Figure [18] shows examples of continuum fits on mock spectra. At z = 2, the estimated continuum is practically 
indistinguishable from the true continuum. The discrepancy between the true and estimated continua increases with 
redshift. At z = 4, the continuum fits underestimate the true continuum by almost 12% on average. 

Let us denote by C es t the estimated continuum level, by Ct ruo the true continuum, and define AC = Ct rU c — C es t- 
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Fig. 18. — Examples of continuum fits to mock quasar spectra. Mock ESI spectra of quasars from redshift z = 2 to z = 4.5 are shown at 
intervals Az = 0.5 successively in the panels. In each case, the thick curve shows the known true continuum level and the blindly estimated 
continuum is indicated by the thinner horizontal curve. At z = 2, the estimated continuum is practically indistinguishable from the true 
continuum. The discrepancy between the true and estimated continua increases with redshift, as quantified in Figure [7] At z = 4, the 
continuum fits underestimate the true continuum by almost 12% on average. In the z = 4.5 panel, we show corresponding spectra, but for 
the background hydrogen photoionization rate, T bkg , half the fiducial value (blue) and twice the fiducial value (red). 

In Figure [71 we plot (AC/Ct rU c) vs. z, where the average is calculated from our 10 continuum fits to mock spectra 
at each redshift. The la uncertainty in this mean is estimated similarly to the error in our (F) measurement (c.f. 
§3.3p . Specifically, each spectrum was broken in segments of 2000 pixels, resulting in 33 chunks at z = 2 and 42 at 
z = 4.5. For each segment i, we calculated the mean (AC/Ct ruo )i, and for each redshift, the standard deviation of 
these segment means, CjAc/Ct,™);- The error 011 the mean at a given redshift is then given by 



r (AC/C tr , 



r (AC/C true )i 



N 



(Bl) 



where the N is the number of segments. The redshift trend is well approximated by the best-fit power law 

AC/Ctmo = 1-58 x 10~ 5 (1 + z) 5 63 



(B2) 

over the redshift range covered. We use this fit to correct our T c g measurement for continuum errors in i)3.4l 

Let the superscripts true and est refer to true values and values estimated from the continuum fitting procedure, 
respectively. Then 
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where the second-to-last step holds because the local transmission can be assumed statistically independent of the 
continuum correction. Combining equations IB31 IB41 and IB51 we find 
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is negligible, as can be seen from Figure [71 and so the statistical error on r*™ is 
just that on rfg'. We discuss the systematic error on our continuum correction in the next section. 



The statistical error on (AC/C tl 
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Uncertainty on our Continuum Correction 

The continuum correction obtained above depends on the properties of the mock spectra we generated. In particular, 
it is a function of the assumed T bkg . The z — 4.5 panel of Figure [T8l for example, shows the mock spectra obtained by 
setting T bkg equal to half and twice the fiducial value, for a fixed line of sight through the IGM. We have used values 
of T bkg that were estimated from previous measurements of T e g (Schaye et al. 2003 for the Bolton et al. 2005 points; 
Fan et al. 2002 for the Meiksin & White 2004 point), which were generally less precise than ours and to which no 
continuum corrections were applied. 12 Our estimated continuum correction may therefore not be self-consistent. 

To estimate by how much our continuum bias correction may be biased, fix some redshift and let Tg Cak be the 
characteristic optical depth of the transmission peaks (which are the points through which we fit a spline to estimate 
the continua; see §3.4j) for the true background photoionization rate F . Let r pcak be the corresponding quantity for 
some other photoionization rate T. Neglecting redshift-space distortions, r poak = Tg eak (r°/T) (e.g., Eq. |2"3"]) . If we use 
an arbitrary value T in generating our mock spectra, then the continuum level estimated from the mock spectra, C os t, 
will differ from the correct value that would have been obtained had we assumed the true r , C cst o : 



C es t \ F c 



Cest.O / 



= e —T^ = exp hrP oak (r /r - 1)] « 1 - r P cak (r /r - 1), (B7) 

peak e T o 



where the last step is valid if the argument of the exponential is <gc 1 . Define Ar c ff = r*g uo (ro) — r*g uo (r), where 
the argument of r e ff indicates the background photoionization rate that is assumed in correcting for continuum bias. 
Substituting Equation (|B6|) . simplifying, and using Equation (|B7[) then gives 

(C es t/Cest,o)(r) 



Ar off 

— r^r = In 



. (C es t/Cest,o)(ro) 



true _ ln[l -r pcak (r/r - 1)] r poak 

r eff - Zt^Tc TtrTIe-i 1 " 1 /^)- (B«) 
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An idea of the magnitude of this fractional error is obtained by assuming that the continuum correction we have 
calculated using values of T bkg from the literature is close to exact. Since the estimated continua pass through 
the transmission peaks, (AC/C tIue ) ~ T pcak . Using the fits in Equations (JUJl and (|B2)) for t* t s uc and (AC/C tIue ), 
respectively, 



At, 



—true 

T eff 



0.0088(1 + z) 
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(B9) 



iBolton et all (|200l performed a detailed analysis of the sources of uncertainty in their determination of T bkg from 
the flux decrement method, concluding that it is known within « 50% at 2 < z < 4, giving 

< 0.0044(1 + Z) 1 - 71 . (B10) 

T eff 

Thus, our continuum-corrected t c s measurement may be off by as much as « (3,5,7)% at z = (2, 3, 4). 13 Note, 
however, that the continuum correction is expected to be a smooth function of redshift, so that uncertainties on it are 
unlikely to affect any conclusion that we may draw about the presence of a narrow feature in the redshift evolution of 

T c fi ■ 

The continuum correction (AC/Ct rU e) ~ TQ Cak also depends on the cosmological parameters through the expression 
for the optical depth, 

12 



„2 



fhya 1 R{T)n ml n e ^ fl 2 b 



where the variables are defined as in equation IB1 11 and we have neglected the cosmological constant a t the redshifts 
z > 2 of interest. However, these cosmological parameters are known much more precisely than T bkg (e.g. JSpergel et all 

l2007h . We see from equation IB 111 that To eak also depends on the temperature T of the IGM, but this dependence is 
degenerate with that on T bkg . The uncertainties on F bkg quantified bv lBolton et al. (2005) include a contribution from 
the uncertainty on T . The degeneracy therefore implies that the error estimates above on our continuum-corrected 
measurement are actually overestimates. 

A more accurate and self-consistent determination of the continuum correction can in principle be achieved by 
iteratively re-estimating it, at each step inferring T bkg from the data using the current correction, and using it for the 
next iteration until convergence is attained. We defer this procedure to future work, in which we will accurately infer 
Y bkg . We here simply emphasize the uncertainty, estimated above, on continuum correction. 
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